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FOREWORD 


At  the  present  time,  thousands  and  thousands  of  people  around  the  world  deal  with  ice, 
snow  and  permafrost.  They  are  scientists,  educators,  engineers,  navigators,  meteorologists 
and  others.  While  a  small  fraction  of  these  people  contribute  to  the  knowledge  base  in  ice 
physics,  all  of  them  use  knowledge  from  it  frequently.  Moreover,  successful  anplied 
research  is  based  upon  fundamental  science — one  more  reason  for  ice  specialists  to  have  a 
textbook  on  ice  physics  on  their  desks. 

The  first  modern  ice  physics  text  was  Fletcher's  book  on  The  Chemical  Physics  of  Ice  ( 1 970). 
Fletcher's  book  is  in  typical  textbook  format:  it  is  reasonably  brief  and  easy  to  understand. 
He  touched  on  a  few  of  the  most  important  topics,  but  not  all  of  them. 

The  most  comprehensive  book  on  ice  physics  to  date  was  published  by  Hobbs  in  1974. 
Hobbs  considered  almost  all  of  the  basic  aspects  of  ice  as  understood  at  that  time.  Moreover, 
he  described  and  compared  several  (sometimes  opposing)  viewpoints.  This  fundamental 
and  rather  large  (837  pages)  book  is  commonly  known  as  the  "Ice  Bible"  by  specialists  in  the 
field.  In  1974  and  1975,  two  CRREL  Monographs  on  ice  were  produced  by  John  Glen.  These 
were  briefly  and  clearly  written  and  reviewed  almost  all  ice-related  subjects.  This  work  was 
(and  in  some  respects  still  is)  a  magnificent  introduction  to  ice. 

Finally,  in  198 1  Maeno  wrote  a  simple,  popular  hook  for  the  express  purpose  of  attracting 
people's  attention  to  the  subject. 

During  the  past  20  years,  a  significant  amount  of  new  experimental  and  theoretical  work 
has  appeared,  dramatically  changing  our  v  iews  on  ice  physics.  As  a  result,  we  are  now  able 
to  formulate  physical  laws  using  more  simple  and  direct  methods.  We  have  found  some  of 
the  physical  models  used  in  the  past  to  be  completely  wrong.  The  physics  of  ice  is  a  much 
better  developed  subject  than  it  was  20  years  ago. 

For  the  above  reasons,  1  feel  the  time  is  ripe  for  a  contemporary  book  on  ice  physics, 
incorporating  the  known  and  proven  with  almost  20  years  worth  of  material  not  covered 
by  previous  works. 

I  have  tried  to  prepare  a  "readable"  book,  and  not  one  that  requires  the  reader  to  be  a 
uniquely  educated  person.  It  is  my  intent  to  present  the  material  in  such  a  way  that  any 
reader  attracted  by  tire  title  lee  Physics  will  be  able  to  comprehend  it.  This  is  quite  difficult 
for  a  book  dedicated,  not  to  a  particular  field  of  knowledge,  but  to  a  specific  material. 
Indeed,  for  >ce  it  means  we  have  to  consider  a  wide  variety  of  subjects,  including  quantum 
chemistry,  solid  state  physics,  the  theory  of  elasticity,  ionic  conductivity,  synchrotron  X-ray 
topography,  crystal  growth,  the  physics  of  surfaces  and  more. 

lire  primary  goal  is  to  produce  as  simple  a  book  as  possible  without  sacrificing  scientific 
accuracy.  Experimental  facts,  physical  ideas  and  tins  les  will  be  strongly  organized  and 
bound  together  cohesively.  Tire  reader  will  be  in  tun  ced  to  a  wide  variety  of  material  on 
a  step  by  step  basis.  Then  the  picture  will  be  whole. 

To  accelerate  publication,  this  book  will  appear  first  in  the  form  of  a  series  of  joint 
CRRFL-Dartrnouth  reports,  later  to  he  published  in  CRREL's  Monograph  series,  on: 

1 .  The  structure  of  ordinary  ice 

Fart  1:  "Ideal'  structure  of  ice  Ice  crystal  lattice 
Fart  II:  Defects  in  Ice 
Volume  1:  Foint  defects 
Volume  2:  Dislocations  and  planar  defects 

2.  Electrical  properties  of  ice 

Fart  I:  Conductivity  and  dielectric  permittivity  of  ice 
Fart  II:  Advanced  topics  and  new  physical  phenomena 
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3.  Optical  properties 

4.  Electro— optical  effects  in  ice 

5.  Thermal  properties 

6.  Mechanical  properties  of  ice.  Elasticity  and  anelastic  relaxation.  Plastic  proper¬ 

ties.  Fracture  of  ice 

7.  Electromechanical  effects  in  ice 

8.  The  Surface  of  ice 

Part  I:  Experimental  results  on  the  structure  and  properties  of  tire  ice  surface 
Part  II:  Theoretical  models  of  the  ice  surface 

9.  Other  forms  of  ice  and  their  properties 

10.  Ice  in  space 

11.  Ice  research  laboratories 

The  reports  will  be  prepared  in  a  sequence  convenient  to  the  author.  The  present 
report  is  the  fifth  in  the  series. 


NOMENCLATURE 

electric  potential 

a 

length  of  hop 

<Ps 

surface  electric  potential 

-1 
®  1 

-i 

JC  ;; 

length  of  screening  by  Bjerrum  defects 

k 

absorption  coefficient 

length  of  screening  by  ions 

Boltzman  constant 

-l 

U_‘ 

screening  length  when  1  e,<ps  1  » k^l 

L 

thickness  of  the  surface  layer 

a0 

fraction  of  surface  molecules  in  the  liquid¬ 

surface  charge  density 

like  layer  with  dipole  moments  (i.e.,  pro¬ 
tons)  oriented  towards  the  surface 

^sa 

surface  charge  absorbed  on  the  ice  inter¬ 
face 

B 

magnetic  induction 

tit 

dipole  momentum  of  a  water  molecule 

Cs 

Cm 

concentration  of  adsorbed  molecules  on 
the  surface 

maximum  concentration  of  adsorbed 

n 

nB 

refractive  index 

concentration  of  Bjerrum  defects 

molecules 

»H20 

water  molecule  concentration  in  ice 

d 

thickness  of  tire  quasi-liquid  layer 

charge  carrier  concentration,  i  =  1,2, 3,4 

D; 

diffusion  roefficient 

i 

refractive  coefficient  of  ice 

angles 

«w 

refractive  coefficient  of  water 

£a 

heat  of  adsorption 

^ion 

concentration  of  ions 

E,  E 

electric  field  strength 

V 

frequency  of  ion's  hops 

<-7 

electric  charge  of  i'1'  type  of  charge  carrier. 

P 

pressure 

i  =  1, 2,3,4 

P 

polarization  vector 

to 

dielectric  permittivity  of  vacuum 

Pc 

—  *• 

electronic  polarization  vector 

es 

static  d’electric  permittivity  of  ice,  Es  =  100 

Pm 

molecular  polari;  ation  vector 

(7  =  -10°C) 

Q 

electric  charge 

£s 

activation  energy  of  surface  conductivity 

P 

rl|/rx 

E' 

effective  electric  field  acting  on  a  molecule 

pW 

density  of  the  space  charge 

high-frequency  dielectric  permittivity  of 

Pl 

water  density 

ice,  e„  ~  3.2 

r± 

reflection  coefficient  for  light  polarized 

<I> 

=  3.85  /cuTr00 

perpendicular  to  the  plane  of  incidence 

<l> 

wetting  angle 

r\\ 

reflection  coefficient  for  light  polarized 

G 

Gibbs  free  energy 

parallel  to  the  plane  of  incidence 

,sr 

gyroinagnetic  constant 

r0 

constant 

Y 

surface  energy 

s 

surface  area 

Ysl 

surface  energy  of  ice/water  interface 

On 

Bjerrum  defect  conductivity 

Yvi. 

surface  energy  of  vapor/ water  interface 

Oil 

bulk  conductivity 

Tvs 

surface  energy  of  va.  >or/ice  interface 

Gor 

ionic  conductivity 

/ 

current 

a„ 

low-frequency  bulk  conductivity 

n 

viscosity 

o.. 

surface  conductivity,  12  ] 

<rss  specific  bulk  conductivity  of  the  ice  sur¬ 
face  layer,  ll'rrr1 
5S|  adhesive  strength 

S0  entropy  of  fusion  per  unit  mass 

high-frequency  bulk  conductivity 
T  temperature 

t  time 

x  correlation  time  of  molecular  motion 

mean  lifetime  of  a  molecule  adsorbed  on 
tire  surface 

tu  Debye  relaxation  time 

Tjjg  Debye  relaxation  time  of  tire  surface  layer 
of  ice 


tms  Maxwell  relaxation  time 

Tm  melting  point  of  ice 

n  velocity 

.v  coordinate 

co  angular  frequency 

to,,  resonance  frequency 

Wj  electron  work  function  of  ice 

Wm  electron  work  function  of  a  metal 

p  chemical  potential 

p,  charge  carrier  mobility,  /  =  1, 2,3,4  forH30+, 

OH-  ions,  D-  and  L-defects  respectively 


The  Surface  of  Ice 


VICTOR  E.  PETRENKO 

INTRODUCTION 


Since  childhood,  we  have  accustomed  our¬ 
selves  to  matter  being  either  solid  or  liquid  or  gas. 
Later,  we  become  acquainted  with  the  fourth  kind 
of  matter  called  plasma,  which  is  very  rarefied  and 
has  an  extremely  large  electrical  conductivity. 

Can  you  imagine  a  simultaneous  combination 
of  certain  properties  of  solid,  liquid  and  plasma 
matter?  It  is  present  on  tire  surface  of  the  ice  upon 
which  we  walk  each  winter,  unsuspecting  that  it 
possesses  such  unique  characteristics. 

When  we  bump  into  a  sharp  corner  or  hurt 
ourselves  with  a  needle  or  cut  a  finger  with  a  knife, 
we  confirm  that  surfaces  of  solids  are  indeed  solid. 
This  is  true  in  so  many  cases  that  it  necessarily 
becomes  a  rule. 

In  tire  case  of  ice,  however,  some  simple  obser¬ 
vations  make  us  doubt  that  it  has  a  rigid  and  per¬ 
manent  surface.  For  example,  if  we  bring  two  ice 
balls  together,  then  right  before  our  eyes,  a  small 
"neck"  starts  to  grow  between  them  (Fig.  1).  In 
time,  further  growth  will  lead  the  two  balls  to 
become  one.  hr  this  simple  experiment,  some  sub¬ 
stance  can  be  considered  to  flow  between  tne  balls 
through  lire  "neck."  Since  ice  in  bulk  form  is  ap¬ 
parently  solid  (have  you  ever  fallen  on  ice?),  it  is 
most  likely  the  surface  that  flows. 


a.  On  contact.  b.  319  minutes  later. 

Figure  1.  Growth  of  a  neck  between  two  ice  spheres  0.5 
mm  in  diameter  at  ~6’C  (from  Hobbs  1974). 


In  a  process  similar  to  the  ice  balls'  unification, 
airy,  fluffy  snowflakes,  after  lying  together  for 
some  time,  suiter  first  into  dense,  heavy  snow,  and 
next  into  solid  ice.  It  is  this  process,  going  on  for 
thousands  and  thousands  of  years,  that  has  formed 
the  kilometer-thick  ice  sheets  in  Antarctica  and 
Greenland  (mass  transport  through  the  vapor 
phase  and  a  plastic  flow  are  also  involved  in  this 
sintering  process). 

Evidently,  the  first  documented  guess  for  the 
existence  of  a  special  layer  on  the  surface  of  ice  was 
the  work  by  Faraday  (1859).  On  the  grounds  of  his 
long-standing  experiments  on  the  adhesion  of  ice, 
Faraday  called  this  layer  a  "liquid-like"  layer.  An¬ 
other  common  expression  is  "quasi-liquid"  layer. 
Now,  the  existence  of  some  special  layer  on  the  ice 
surface  has  been  proven  by  many  experimental 
techniques.  Hundreds  of  papers  are  dedicated  to 
the  investigation  of  the  structure  of  tills  layer  and 
its  physical  properties.  The  most  essential  contri¬ 
butions  to  the  store  of  our  knowledge  about  this 
liquid-like  layer  were  made  during  the  past  two 
decades. 

In  the  present  repoit,  I  shall  describe  the  results 
of  the  experimental  research  that  provides  direct 
information  on  die  structure  and  thickness  of  this 
layer:  X-ray  diffraction,  proton  channeling,  ellips- 
ometry  and  NMR.  Then,  numerous  results  of  stud¬ 
ies  of  the  physical  properties  of  tli-a  ice  surface  will 
be  considered.  An  analysis  of  these  properties  by 
means  of  comparing  various  models  also  allows  us 
to  obtain  indirect  information  about  the  structure 
and  thickness  of  the  special  layer  on  the  ice  surface. 
We  shall  not  consider  ice  friction  here,  since  sever¬ 
al  good  reviews  on  this  topic  have  been  published 
recently  (see,  for  instance,  Aharon  et  al.  [1988]  and 
Colbeck  [1992]).  Nevertheless  when  all  the  reports 
are  united  into  a  book,  a  chapter  on  ice  friction  will 
be  included. 

In  conclusion,  major  known  theoretical  models, 
describing  the  structure  and  properties  of  the  liq¬ 
uid-like  layer,  will  be  discussed. 


PART  I.  EXPERIMENTAL  RESULTS  ON  THE  STRUCTURE 
AND  PROPERTIES  OF  THE  ICE  SURFACE 


X-RAY  DIFFRACTION 

When  considering  the  structure  of  the  quasi¬ 
liquid  film  on  the  ice  surface,  we  must  first  ans¬ 
wer  the  question:  is  the  structure  of  the  film  or¬ 
dered  or  not?  That  is,  does  the  structure  have  a 
long  range  order?  And  if  so,  then  is  this  structure 
different  from  that  of  ice?  It  seems  that  the  easiest 
way  to  answer  this  question  is  to  use  any  known 
X-ray  diffraction  technique.  However,  owing  to 
lire  extremely  small  thickness  of  the  film,  it  is  dif¬ 
ficult  to  differentiate  its  contribution  to  the  dif¬ 
fraction  from  that  of  the  ice  bulk.  I  am  aware  of 
only  one  such  attempt,  where  Kouchi  et  al.  (1987) 
investigated  X-ray  diffraction  at  the  surfaces  of 
polycrystalline  and  monocrystalline  ice  in  a  tem¬ 
perature  range  of  -0.5  to  -10°C.  Glancing  beams 
were  used  (angle  of  incidence  =  2°)  to  make  the 
contribution  of  the  surface  as  large  as  possible. 
Tire  Ka  line  of  the  Cu  X-ray  spectrum  was  used  in 
the  experiments. 

Figure  2  shows  X-ray  diffraction  patterns  of 
tire  surface  of  the  polycrystalUne  ice  sample  ob¬ 
tained  at  -9.3,  -2.1  and  -LHC'.  At  T  =  -9.3  and 
-2.1°C,  the  diffraction  peaks  that  are  characteris¬ 
tic  for  the  crystal  are  clearly  visible,  although 
their  magnitude  decreases  at  -2.1°C.  At  T  =  -1.1 
°C.  those  peaks  virtually  vanish,  being  re¬ 
placed  by  a  blurred,  circular  "halo"  cen¬ 
tered  at  20  =  26°.  Tlr.s  is  a  clear  evidence 
that,  at  this  temperature,  the  long-range 
order  in  the  arrangement  of  water  mole¬ 
cules  in  the  surface  layer  vanishes. 

Since  the  intensity  of  the  "halo"  can  be 
used  to  estimate  the  disordered  layer's 
thickness,  it  is  interesting  to  follow  the 
temperature  dependence  of  this  intensity, 
which  is  illustrated  in  Figure  3.  First,  we 
can  conclude  from  this  figure  that  the  lay¬ 
er's  thickness  is  independent  of  the 
crystallographic  orientation  of  the  sur¬ 
face.  Second,  the  layer  ceases  to  be  detect¬ 
able  by  this  method  at  T  <  -2°C.  Unfortu¬ 
nately,  Kouchi  et  al.  did  not  attempt  to 
estimate  an  absolute  magnitude  of  the 
layer's  thickness. 


Figure  2.  Intensity  of  the  beam  diffracted  from  the 
polycrystalUne  ice  surface  as  a  function  of  the  angle 
between  incident  and  reflected  beams  (after  Kouchi  etal. 
1987). 


temperature  (°C) 

Figure  3.  Changes  in  intensity  of  diffraction  halo  at  2x)  =  28°  with 
temperature  (after  Kouchi  et  al.  1987). 


PROTON  CHANNELING 


Another  interesting  experimental  method  of 
determining  the  thickness  and  structure  of  the 
quasi-liquid  layer  was  used  by  Golecki  and  Jac- 
cara  (1977, 1978).  This  is  the  so-called  channeling 
effect,  discovered  by  Oen  and  Robinson  (1962). 
The  channeling  effect  manifests  itself  in  the  fact 
that,  if  a  beam  of  heavy  charged  particles  (usually 
protons  or  a  pa ’•tides)  travels  parallel  to  some 
crystallographic  direction,  then  the  particles  of 
that  beam  experiencing  low-angle  scattering  pen¬ 
etrate  the  crystal  at  very  long  distances  and  can 
stay  inside  the  channels  formed  by  atomic  rows. 
Tire  magnitude  of  this  effect  is  usually  expressed 
in  terms  of  the  ratio  of  the  flux  of  backscattered 
particles  when  tire  direction  of  the  flux  coincides 
with  one  of  the  crystallographic  directions  to  that 
when  the  direction  is  arbitrary. 

The  backscattering  can  take  place  only  at  nu¬ 
clei  of  atoms  that  are  much  heavier  than  the  parti  ¬ 
cles  in  the  beam  and  resembles  the  classic  Ruther¬ 
ford  scattering.  The  diffraction  of  such  particles  at 
the  crystal  lattice  is  insubstantial,  since  the 
characteristic  de  Broglie  wavelength  is  several 
orders  of  magnitude  less  than  the  lattice  constant. 

Most  of  tire  backscattered  particles  experience 
only  a  single  strong  interaction  with  a  lattice 
atom.  Using  the  theory  of  Rutherford's  scatter¬ 
ing,  we  can  calculate  the  energy  of  the  reflected 
particle  as  a  function  of  the  scattering  angle.  De¬ 
partures  from  this  predicted  value  contain  infor¬ 
mation  about  the  depth  at  which  the  scattering 
took  place  (more  exactly,  about  the  path  the  parti¬ 
cle  has  traveled  back  arid  forth  within  tire  matter). 
The  loss  of  the  kinetic  energy  of  the  particles  is 
attributable  to  their  interaction  with  the  valence 
electrons  and  is  almost  independent  of  the  direc¬ 
tion  of  the  motion.  Having  "calibrated"  such  loss¬ 
es  for  a  particular  material  (or  having  calculated 
them  theoretically),  we  can  establish  the  depth  at 
which  tire  scattering  took  place  with  a  rather  high 
accuracy. 

The  channeling  is  extremely  sensitive  to  the 
degree  of  perfection  of  a  crystal  structure  and 
therefore  can  be  efficiently  used  for  measure¬ 
ments  of  the  magnitude  of  thermal  vibrations,  the 
degree  of  disordering,  the  concentration  of  point 
defects  and  so  forth.  One  of  the  comprehensive 
reviews  on  the  subject  is  that  by  Gemmel  (1974). 

In  1977  and  1978  this  effect  was  explored  by 
Golecki  and  Jaccard  for  the  investigation  of  the 
(0001)  surface  of  ice  crystals.  The  authors  mea¬ 
sured  the  backscattering  (angle  150°)  of  100-keV 
promns  from  the  surface  of  monocrystalline  ice  in 


a  wide  temperature  range  from  -2  to  -i30°C. 
They  found  that,  at  temperatures  higher  than 
-40°C,  the  surface  layer  is  characterized  by  much 
higher  disordering  in  relation  to  tne  bulk.  Analy¬ 
sis  of  the  behavior  of  the  disordering  process 
prompted  Golecki  and  Jaccard  to  conclude  that 
the  disorderng  is  caused  by  a  significantly  high¬ 
er  magnitude  of  thermal  vibrations  of  oxygen 
adorns  near  the  surra/  i  ather  than  the  absence  of 
symmetry  in  the  structure.  Thus,  at  T  -  -1.8°C  the 
average  amplitude  of  the  thermal  vibrations  of 
oxygen  atoms  near  the  surface  exceeds  that  of  tire 
bulk  by  3.3  tin.es.  Because  of  that,  Golecki  and 
Jaccard  suggest  not  using  the  term  'liquid-like' 
layer  at  all. 

To  estimate  the  thickness  of  the  special  surface 
layer  of  ice,  Golecki  and  Jaccard  calculated  an 
equivalent  thickness  of  an  amorphous  layer  (to¬ 
tally  disordered)  that  would  result  in  the  same 
magnitude  of  scattering.  The  results  of  such  cal¬ 
culations  are  depicted  in  Figure  4.  It  is  worth 
noting  that  the  half-thickness  of  the  layer  f/2  ob¬ 
tained  is  extremely  large.  The  temperature  de¬ 
pendence  of  the  thickness  t  was  found  to  be 

<f(nm)  =  (94  ±  17)  -  (54  ±  14)  log(273  -  T).  (1) 


273  -T (K) 

Figure  4.  Equivalent  "amorphous''  thickness  of  the 
disordeied  region  at  the  ice  surface,  deduced  from  the 
measurements  with  the  aid  of  Meyer's  (1971)  theory. 
The  solid  hue  corresponds  to  least-squares  fit.  The 
density  of  the  surface  region  was  assumed  equal  to  that  of 
ice  (0.917 g  cm  d)  (after  (Joic'ki  and  jaccard  1978). 
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OPTICAL  ELLIPSOMETRY 

Ellipsometry  is  an  optical  reflectance  tech¬ 
nique  in  which  the  change  of  the  state  of  polariza¬ 
tion  of  light  upon  reflection  from  a  surface  is 
measured.  Even  though  the  method  dates  back  to 
Drude  (1889),  ellipsometry  has  recently  attracted 
considerable  attention  because  of  a  wide  variety 
of  applications  in  modem  solid  state  optics,  espe¬ 
cially  in  tire  physics  of  surfaces,  interfaces  and 
thin  films.  In  the  study  of  the  ice  surface,  optical 
ellipsometry  was  used  first  by  Beaglehole  and 
Nason  (1980),  and  later  by  Furukawa  et  al. 
(1987a, b). 

The  change  of  the  state  of  polarization  upon 
reflection  can  be  expressed  in  terms  of  the  ratio  of 
two  complex  reflection  coefficients  p|  and  rj  for 
light  polarized  parallel  and  perpendicular  to  the 
plane  of  incidence,  r  being  the  ratio  of  refected 
and  incident  electric  field  strength.  The  complex 
quantity 

p  =  — •  =  tan\|(  exp  (i'A)  (2) 

defines  the  two  ellipsometric  angles  A  and  \y.  The 
two  angles  completely  determine  tire  two  optical 
constants  ti  (refractive  index)  and  k  (absorption 
coefficient)  of  an  isotropic  reflecting  medium. 
Mathematical  dependence  of  n  and  k  on  A  and  vp, 
calculated  from  Fresnel's  formula,  was  found,  for 
example,  by  Bootsma  and  Meyer  (1969). 

The  application  of  ellipsometry  in  surface 
physics  is  favored  because  in  suitable  cases  a  lay¬ 
er  with  average  thickness  as  small  as  Vioo  of  a 
monolayer  can  be  detected,  if  in  tire  layer  u  and  k 
are  different  from  that  of  the  bulk. 

Ellipsometric  data  are  usually  analyzed  with 
computers  because  the  two  measured  variations 
5A  and  8rp  are  not  sufficient  to  determine  all  pa¬ 
rameters  (optical  constants,  thickness)  even  for  a 
single  surface  film. 

Beaglehole  and  Nason  (1980)  employed  ice 
monocrystals  grown  following  the  Bridgeman 
method  using  twice-distilled  water.  The  ice  sur¬ 
face  with  the  orientation  (0001)  (basal  plane)  or 
(1010)  (prismatic  plane)  was  milled  with  a  warm 
razor  immediately  before  ellipsometrical  tests. 
Water  forming  simultaneously  on  the  ice  surface 
was  removed  by  the  same  razor,  resulting  in  a 
surface  of  unruffled  cleanness. 

Measurements  in  the  0  to  -30°C  temperature 
range  rev  aied  a  considerable  distinction  in  tire 
behavior  of  tire  surface  layer  at  the  basal  and  pris¬ 
matic  planes.  Assuming  that  the  reh  action  coeffi¬ 
cient  of  the  surface  layer  is  the  same  as  that  of 


water,  Beaglehole  and  Nason  found  that  the  film's 
thickness  on  the  prismatic  surface  ( d )  was  about 
13C  A  at  T  =  -1.5°C,  decreasing  slowly  to  a  mono- 
atomic  thickness  at  T  =  -10°C.  They  found  a  sur¬ 
face  layer  on  the  basal  plane  only  at  temperatures 
greater  than  -0.3°C. 

Furakawa  et  al.  (1987a,b)  significantly  ad¬ 
vanced  use  of  the  ellipsometry  technique  for  de¬ 
termining  tire  transition  layer  thickness  on  the  sur¬ 
face  of  ice.  First,  they  used  unproved,  more  perfect 
ice  crystals,  grown  from  twice-distilled  and  de¬ 
gassed  water.  The  degree  of  purity  of  ice  plays  an 
important  role,  since  even  a  small  salt  concen- 
'■ i .  it  ion  can  sizably  alter  the  thickness  of  the  transi¬ 
tion  (quasi-liquid)  layer  on  tire  ice  surface.  Second, 
they  examined  atomically  flat  (0001)  and  (1010) 
surfaces  of  the  natural  facets  of  a  so-called  nega¬ 
tive  crystal.  A  negative  crystal  is  formed  by  pump¬ 
ing  water  vapor  out  of  a  small  void  inside  a  large 
monocrystal.  Besides,  they  made  their  ellipsomet¬ 
ric  measurements  under  conditions  of  equilibrium 
between  ice  and  its  vapor. 

Finally,  Furakawa  et  al.  interpreted  their  data, 
not  using  tire  assumption  that  the  refraction  coeffi¬ 
cient  of  the  surface  film  (>;,)  is  the  same  as  that  of 
water  (riw),  but  calculating  it  as  a  fitting  parameter. 
They  found  that  n1  ~  1 .330  for  both  faces  (0001 )  and 
(1010),  which  is  closer  to  water,  nw  =  1.3327,  than 
to  ice,  n,  =  1.3079.  This  implies  that  the  density  of 
the  layer  is  quite  close  to  that  of  water. 

Figure  5  shows  temperature  dependence  of  the 
transition  layer  thickness  as  found  by  rurakawa  et 
al.  from  their  ellipsometrical  measurements.  As 
follows  from  this  plot,  the  thickness  of  the  transi¬ 
tion  layer  drastically  rises  as  temperature  increas- 
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Figure  5.  Temperature  dependence 
of  the  relative  phase  change  A  and 
the  thickness  of  liauid-like  layer  d 
( after  Furukawa  et  al.  1987a). 
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cs  above  I  l  rximlmt  lr'H  nm  A  dilferenn- 
between  tbt'  film  s  dm  kin->M  ■>  on  the  (iMH>l  )  t.u  <• 
and  the  ( i 0 10)  (.lie  is  espix  i.illv  notu  cable  in  the 
temperature  i.mge  (rom  1  to  1  G  Optical  e!lip 
sometrv  can't  detei  t  the  pieseiu  e  ot  a  liquid  ulm 
on  an  a  e  surtax e  at  temperatures  below  -1  ( 

Wrv  reeenth  ,  Elbaum  et  al  ( l1*4"))  performed 
optical  measurements  ot  a  liquid  like  laser  using 
a  .simplified  elhpsometrv  technique.  which  was 
termed  Brewster  reflectometrv  lite  idea  was  to 
measure  the  refr.u  live  index  profile  near  a  plane 
interface  between  uniform  media.  It  a  collimated 
light  l*eam  wi'.'i  a  vector  /.  lving  in  the  plane  ot 
incidence  falls  at  the  Brewster’s  angle,  then  the 
intensity  of  the  reflected  light  will  be  zero  for  an 
interface  between  two  homogeneous  media.  Bui, 
since  the  interface  has  a  finite  thickness,  the  re¬ 
flectivity  is  slightly  different  from  zero.  Elhaum 
et  al.  assumed  that  the  interface  is  composed  of 
water  vap  xr  or  its  mixture  with  air  having  the 
refraction  coefficient  n  =  1,  a  liquid  film  of  thick¬ 
ness  L  and  n  =1.333  (same  as  water)  and,  finally, 
ice  with  n  -  1.309.  Given  such  conditions,  the 
problem  of  the  intensity  of  the  reflected  light  can 
bo  solved  exactly,  and  they  used  this  solution  to 
calculate  the  thickness  of  the  liquid  film  on  the  ice 
surface  from  the  measurements  of  intensity  of  a 
polarized  light  reflected  at  the  Brewster's  angle. 
Though  this  method  has  a  lower  sensitivity  com¬ 
pared  to  conventional  ollipsometry  (L  2  12  A, 
Elbaum  et  al.  [1993]),  the  described  work  was  the 
first  where  the  effect  of  air  on  the  thickness  and 
conditions  of  existence  of  the  quasi-liquid  layer 


Figure  6.  Collected  results  on  the  liquid-like  film  thick¬ 
ness  from  four  reflection  experiments.  The  inset  shows 
the  same  results  in  smaller  thickness  range  (after  Elbaum 
et  ai  l  993). 


u.w  ohseivcd  and  mv vMig.itod  In  the  aheeme  of 
air.  the  quaM-liquid  him  could  not  be  detec  ted  at 
the  prismatic  tan-  ol  ue  within  the  experimental 
erior  (12  A)  At  the  came  time,  a  \cr\  thick  him 
was  obseiwd  at  the  basal  lace  (IKHM)  order  the 
same  conditions  (1  ig  <•>)  1  he  presein  e  ot  an  Mgtul 
uamlv  increased  the  liquid  him  thickness  at  the 
<(XKH)  lace 

NUCLEAR  MAGNETIC  RESONANCE 

Clifford  (1907).  using  Nuclear  Magnetic  Keso 
nance-  (NMR)  measurements  on  ice.  observed  a 
verv  narrow  line  (  0.2  G)  ihat  he  attributed  to 
some  form  of  liquid  water  remaining  in  ice  down 
to  temperature  ot  at  least  12"C.  Kvlivid/.c  et  al. 
(1970)  demonstrated  that  this  narrow  resonant 
line  is  distinct  from  the  NMR  lines  corresponding 
to  protons  in  water  and  bulk  ice  and  is  caused  by 
the  liquid-like  layer  on  the  ice  surface.  Afterwards, 
NMR  was  used  for  studying  the  properties  of  the 
liquid-like  layer  on  the  ice  surface  by  Bell  et  al. 
(1971),  Kvlividze  et  al.  (1974),  Barer  et  al.  (1977) 
and  Ocampo  and  Klinger  ( 1983a, b).  The  most 
comprehensive  information  about  the  physical 
properties  of  a  liquid  like  layer  was  obtained  by 
Mizuno  and  Hanafuza  ( 1987),  to  whom  u  e  are  go¬ 
ing  to  refer  most  often  for  this  reason. 

The  phenomenon  of  NMR  is  an  absorption  of 
electromagnetic  field  quanta  by  nuclei  that  pos¬ 
sess  a  magnetic  moment  and  are  placed  in  an  ex¬ 
ternal  magnetic  field  II  at  resonance  frequency  (0o 

(0 „  =  gll  (3) 

where  g  is  the  giromagnetic  constant. 

The  position  of  the  resonant  line,  its  shape  and 
half-width,  as  well  as  the  spin-lattice  relaxation 
time,  appear  to  be  quite  sensitive  to  the  properties 
of  a  medium  that  surrounds  the  nuclear  spins. 
Namely,  those  properties  are  the  average  distance 
and  relative  position  of  the  spins  in  space,  angular 
velocity  and  diffusion  coefficient  of  molecules  that 
contain  the  resonant  nuclei. 

Readers  who  are  interested  in  this  method 
should  see  the  marvelous  book  by  Abragam  (1961) 
and  also  the  paper  by  Burnett  and  Harmon  (1972). 

The  NMR  technique  applied  to  ice,  water  or  the 
quasi-liquid  layer  allows  us  first  to  distinguish 
three  resonant  lines  associated  with  absorption  by 
hydiogen  nuclei  in  those  three  media.  From  the 
intensity  of  the  liquid-like  layer  line,  we  can  judge 
the  approximate  thickness  of  a  particular  layer. 
Besides,  it  is  possible  to  determine  the  c  orrelation 
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time  ol  rotational  motion  ot  water  molecule-.  and 
a  solt-dittuMim  inefficient  (Mizuno  and  Manatu 
/a  14K7) 

Since  under  normal  londitinns.  the  relative 
portion  of  molecules  in  the  liquid  like  lay  nr  is 
negligibly  small,  it  is  necessary  first  to  increase 
the  "surfaco-to-volume"  ratio  of  tho  ice  in  ques- 
turn  to  observe  NMK  from  this  laver.  IT»e  easiest 
way  to  achieve  this  is  to  pulverize  water  at  a  low 
temperature.  As  a  result,  microscopic  ice  bails  are 
formed  with  a  diameter  about  100  pin  and  a  net 
surface  of  about  3-SxU}2  cm/g.  This  method  was 
used  fo;  instance  hv  Mizuno  and  Hanafuza 
(1487).  Kvlividze  et  al.  (1474)  sjturatc'd  water 
with  air  prior  to  freezing.  This  leads  to  an  ap¬ 
pearance  of  a  great  number  of  microscopic  air 
bubbles  inside  the  ice,  with  a  net  surface  of  the 
order  of  ~  104  cm2/g.  Frost  forming  on  walls  of  a 
Dewar  vessel  filled  with  liquid  nitrogen  has  an 
even  larger  specific  area  of  1  -2x10’’  em2/g  {Kvliv¬ 
idze  etal.  1  ‘>70) .  Placing  spherical  Teflon  particles 
(r  «=  350  A)  (Kvlividze  et  al.  1474),  or  Aerosil 
particles  (r  *  100  A)  (barer  et  al.  1477)  into  water 
before  freezing  produc’d  something  comparable 
to  the  "frost"  net  surface  of  ice  for  NMK  analysis. 
We  have  to  take  into  account,  of  course,  that  tile 
structure  and  properties  of  the  liquid-like  layer 
on  lire  ice/ (solid  dielectric)  interface  might  bo  dif¬ 
ferent  from  ice/(water  vapor  +  air)  interface 
properties. 

Figure  7  shows  NMK  spectra  of  the  liquid-like 
layer  of  ice  measured  at  various  temperatures.  To 


Figure  7.  NMli  spectra  of  the  quasi-liquid  layer  of  ice 
and  water  observed  at  99.5  MHz.  The  signal  of  liquid 
water  was  taken  al  5‘C.  Notice  the  difference  in  the  line 
width  between  liquid  water  and  the  quasi-liquid  layer 
(after  Mizuno  and  l  lanafusn  1987). 


figure  8.  Intensity  and  line  width  of  NMR  of  quasi¬ 
liquid  layer  against  inverse  temperature  (after  Mizuno 
and  Hanafusa  1987). 


increase  the  signal-to-noise  ratio,  the  spectra 
were  gathered  200  times  at  each  temperature.  The 
broad  signal  from  crystalline  ice  is  not  seen  with¬ 
in  the  range  of  the  20-kHz  observational  frequen¬ 
cy  (Mizuno  and  Hanafuza  1987). 

The  narrow  signal  was  never  detected  at  any 
temperature  if  bulk  ice  was  used  as  a  sample. 

As  is  clear  from  Figure  7,  the  line  width  and 
the  intensity  vary  with  temperature  and  their 
temperature  dependencies  are  shown  in  Figure  8. 
It  should  be  noted  that  the  line  width  of  the  liq- 
uid-like  layer  spectrum  at  -10°C  is  about  7  times 
that  of  ordinary  water  at  5°C. 

Since  the  temperature  dependencies  shown  in 
Figure  8  were  measured  as  temperature  was  in¬ 
creased  from  -100  to  -5°C,  and  since  at  T  £  -10°C 
a  considerable  sintering  diminishes  the  total  ice 
surface,  the  actual  increase  of  the  NMR  signal  at  T 
£  -10°C  could  be  greater  than  the  one  shown  in 
Figure  8. 

The  temperature  dependence  of  the  correla¬ 
tion  time  for  rotational  motion  x  is  shown  in  Fig¬ 
ure  9.  This  time,  which  could  be  interpreted  as  the 
time  it  takes  a  water  molecule  to  be  reoriented  in 
a  liquid-like  layer  at  T  =  -10C,C,  appears  to  be 
closer  to  that  time  for  water  (=10~12  seconds)  than 
for  ice  (2x10  5  seconds  [Jaccard  1959])  at  the  same 
temperature.  The  activation  energy  of  molecule 
reorientation  in  a  liquid-like  layer,  which  can  be 
determined  from  the  slope  of  the  curve  in  Figure 
9,  is  0,291  eV  and  is  also  quite  different  from  the 
corresponding  value  for  the  bulk  ice  (0.575  eV 
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Figured.  Correlation  time  for  rotational  motion 
r  of  quasi-liquid  layer  against  inverse  tempera¬ 
ture  (after  Mizuno  and  llonafusa  1987). 

[Jaccard  1959]).  Therefore,  comparing  the  rota¬ 
tions  of  a  water  molecule  in  water,  the  liquid-like 
layer  and  ice,  we  -ould  conclude  that,  according 
to  this  property,  the  ilquid-like  layer  is  intermedi¬ 
ate  between  ice  and  water,  but  it  nevertheless 
resembles  water  more  closely. 

However,  if  we  try  to  compare  self-diffusion 
coefficients  in  water,  the  liquid-like  layer  and  ice 
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Figure  10.  Diffusion  coefficient  of  quasi-liquid  layer  against 
inverse  temperature  (after  Mizuno  and  llonafusa  1987). 


(see  big.  10),  we  will  come  to  an  opposite  conclu¬ 
sion,  because  this  coefficient  at  0UC  in  the  liquid¬ 
like  layer  turns  out  to  be  -  3xlO-12  n.2/s,  which  is 
closer  to  ice  (-HH4  m2/s)  than  to  water  (2x10  4 
m2/s).  At  the  same  time  the  activation  energy  of 
self-diffusion  in  the  liquid-like  layer  (0.24  eV)  is 
closer  to  that  of  water  (=0.24  eV)  than  to  ice  (0.65 
eV). 

Barer  et  al.  (1477)  studied  NMR  from  a  liquid- 
like  layer  forming  at  the  interface  between  ice  and 
Aerosil  microspheres,  i.e.,  on  a  frozen  aqueous 
suspension  of  Aerosil.  Like  Mizuno  and  Hanafu 
za,  they  tracked  tire  existence  of  tire  liquid-like 
layer  much  further  down  on  the  temperature  scale 
than  measurements  that  could  be  obtained  from  X- 
rays  or  ellipsometry  (see  Fig.  11).  The  results  of 
their  estimations  of  the  liquid-like  layer  thickness 
at  various  temperatures  are  given  in  Table  1. 

Table  1.  Temperature  dependence  of  the 
thickness  of  a  liquid-like  layer  on  the  ice/ 

Aerosil  interface  (after  Barer  et  al.  1977). 
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In  the  same  work  (Barer  et  al.  1977),  the  viscos¬ 
ity  of  a  liquid-like  layer  on  the  ice/quartz  interface 
was  measured.  At  T  -  ~10"C  the  viscosity  ap¬ 
peared  to  be  i)  =0. 17 1’,  which  exceeds  the  viscosity 
of  supercooled  water  at  the  same  temperature  by 
an  order  of  magnitude. 

N„% 


210  230  lb 0  2/0  K 


Figure  II.  Temperature  dependence  of  the  relative 
concentration  of  the  mobile  aqueous  suspension  of 
Aerosil  (after  Oarer  et  al.  1977). 


ELECTRICAL  CONDUCTIVITY 
OF THE ICF SURFACE 


As  wo  have  seen  in  the  previous  sections,  the 
surface  layer  of  ice  differs  from  the  bulk  in  its 
structure,  optical  properties,  molecule  reorienta¬ 
tion  time  and  self-diffusion  coefficient.  It  is  there¬ 
fore  also  natural  to  expect  the  ice  surface  to  be  dif¬ 
ferent  from  the  bulk  in  its  electrical  properties. 

We  can  name  several  possible  things  that  could 
change  electrical  properties  of  ice  near  the  surface. 
The  first  is  an  increase  (or  decrease)  of  the  charge 
carrier  concentrations  within  the  screening  length, 
assuming  that  there  is  a  potential  difference  cp„ 
between  the  surface  and  the  bulk.  Such  a  differ¬ 
ence  may  arise  because  of  differences  in  the  chem¬ 
ical  potentials  of  charge  carriers  in  the  bulk  and  on 
the  surface  (electrochemical  potential  must  be  uni¬ 
form  everywhere)  or  because  of  adsorption  of  elec¬ 
trically  charged  impurities  that  compensate  the 
space  charge  of  ice  at  the  surface,  The  latter  arises 
merely  from  differences  in  solubilities  of  positive 
and  negative  impurity  ions  in  ice.  Second,  charge 
carrier  mobilities  and  the  mechanisms  of  interac¬ 
tion  between  them  might  be  different  in  the  disor 
dered  surface  or  quasi  liquid  layer. 

In  fact,  the  ice  surlace  has  an  extremely  high  not 
electrical  conductivity  that  was  discovered  back  in 
the  lbSOs  and  has  been  studied  intensely  over  the 
past  30  years,  At  temperatures  roughly  above 
-  LTV  m  pure  ice,  the  d  ■. .  conductivity  ol  the  sur 
lace  exceeds  that  el  the  hulk.  A  measurement 


ing  on  the  surface  from  the  bulk  currents.  Usual¬ 
ly,  this  is  done  by  means  of  so-called  guard  rings 
(Fig.  12).  A  guard  r;ng,  as  can  be  seen  from  this 
figure,  is  separated  from  the  central  electrode  by 
a  thin  insulating  ring.  The  guard  ring  and  the 
central  electrode  have  nearly  the  same  potential, 
which  provides  an  orientation  of  electrical  field 
strength  vector  L  (and  consequently  of  currents) 
parallel  to  the  surface  of  the  cylinder.  In  this  set¬ 
up  all  subsurface  currents  flow  through  the 
guard  ring,  and  their  magnitude  is  proportional 
to  the  voltage  drop  across  the  resistor  Ks,  In  this 
schematic,  Rs  and  R|t  must  be  much  smaller  than 
the  impedance  of  the  specimen  Zs  in  the  whole 
range  of  the  investigated  frequencies. 

A  great  majority  of  researchers  have  used  the 
configuration  shown  in  Figure  12,  although 
sometimes  square  ice  samples  wore  used  instead 
of  cylindrical  ones.  Caranti  and  Illingworth 
( l‘W3a)  used  the  same  idea  (making  parallel  to 
the  surface  and  the  use  of  a  guard  ring),  although 
they  employed  a  semi-spherical  electrode  frozen 
into  ice. 

] award  (1%<»,  1%7)  suggested  that  four-point 
electrodes,  arranged  in  pairs  at  the  opposite  sides 
of  a  thin  ice  specimen,  be  used  to  determine  the 
suifaiv  and  bulk  conductivities.  However,  we 
shall  refrain  from  discussing  this  here  because  of 
serious  eritieism  of  this  meihod  expressed  by 
Turner  and  Slow  (llWh). 


of  the  sui  face  conductivity  is  one  ol  the  sim 
plest  ami  most  accessible  ways  to  study  the 
physical  properties  ol  the  nesurlaee.  in  fact, 
a  simultaneous  measurement  ol  the  absolute 
value  ol  l he  surface  conductivity  os  (with  a 
dimension  of  U  ')  and  ot  the  specilic  volume 
electrical  conductivity  ot  the  surface  layer  o*, 
(with  a  dimension  ot  li  'ill  •)  would  allow 
determination  oi  the  layer’s  thickness  I, 

/.  -  n>  .  (4) 

cc.. 

Measurement  of  charge  carrier  mobilities' 
values  [i,  and  their  concentrations  n,  is  essen¬ 
tial  to  understanding  the  nature  and  the 
structure  of  the  quasi-liquid  layer,  as  well  as 
for  determining  the  activation  energies  of  p, 
and  iir  Researchers  in  this  field  have  investi¬ 
gated  the  electrical  conductivity  and  dielec¬ 
tric  permittivity  of  the  surface  layer,  as  well 
as  the  1  lall  effect. 

The  major  methodical  problem  in  these  ex¬ 
periments  is  the  separation  of  currents  flow¬ 
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Figure  /  2.  Typical  configuration  for 


vvrn  measurements  of  tin  face  and  bulk 


conduct  ivili/. 


Table  2.  Surface  conductivity  of  ice  . 
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Camp  el  al.  (1464) 
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Id 

107 
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0.7K 

pure  polycrvstals 

Caranti  and  l.amfn  (1487) 
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di¬ 

•  10“ 

1.45 

pure  single  crystals 

Maidique  et  al.  (1471) 

-10 
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1.24 

pure  single  crystals 

Kuepp  and  Kiiss  (1464) 

-11 
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(1  4-6.5)  x  10  111 

— 

pure  single  crystals 

laecard  (1467) 
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d.i. 

.  1  5x10  "* 

polycrvstalline  ice  him 

C  iii  Miunvski  (1488) 

•Id 

di. 

-  4x10 

de|x>sit,ii  unlo  a  dielectric 
subsirale 

*  i  aliiilated  tmiii  Chr/anutvski's  data  assuming  that  the  electrode's  length  was  5  mm. 
t  Calculated  trom  Chrianowskt's  data  assuming  that  the  .novation  energy  is  1.2  eV 


liven  a  small  amount  of  impurities  in  ice 
strongly  affects  the  surface  conductivity  o„,  This 
accounts  for  the  considerable  scatter  in  the 
experimental  data  obtained  on  a».  Data  on  low- 
frequency  and  d.c.  surface  conductivities  of  pure 
and  doped  ice  at  -10  "C  are  summarized  in  Table 
2.  As  we  can  see  from  this  table,  a,,  even  for  pure 
monoeryslalline  ice,  assumes  values  in  a  wide 
range,  from  10  11  i2"!  (Maeno  and  Nishitnura 
1  W7fS)  up  to  10  g  U  1  (Maidique  et  al.  1071).  Data 
on  activation  energies  L„  for  the  same  crystals 
vary  from  a  minimum  value  of  1.1*4  to  1.52  eV, 
with  an  average  value  of  about  1.3  eV.  liven  a 
slight  doping  of  ice  leads  to  a  considerable  in¬ 
crease  in  as  and,  as  a  rule,  a  decrease  of  Z.v 
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Typical  frequency  dependencies  of  ice  bulk 
and  surface  conductivities  are  plotted  in  Figure 
13.  This  plot  shows  that  the  frequency  depen¬ 
dence  of  the  bulk  conductivity  ctu  follows  De¬ 
bye's  dispersion,  which  is  well  known  for  ice 
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(5) 


where  o()  is  the  low-frequency  limit  of  the  con¬ 
ductivity  (or  d.c.  conductivity)  determined  by 
minority  charge  carriers.  In  pure  ice  at  -10"C, 
those  are  1 1(0 ‘  and  OH  ions.  o„  is  the  high- 
frequency  limit  of  the  conductivity,  equal  to  the 
sum  of  all  partial  conductivities  of  the  charge 
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u  -  bulk  (an) 
x  -  surface  (a,) 

•  -  surface  with  thick  (rost 
■  -  suifaco  after  annealing. 

Figure  (3.  Variation*  of  the  real  part  of cotnitu  Unity  of 
pure  polycryotalline  ice  at  -1U°C  (after  Caranti  and 
Illingworth  1983). 
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carriers.  Under  the  conditions  considered, 
practically  coincides  with  the  conductivity  of  the 
majority  charge  carriers — Bjerrum  defects,  that  is 
L-  and  D-defects.  co  is  the  angular  frequency  and 

is  the  Debye  relaxation  time.  I  refer  readers 
who  are  more  interested  in  the  electrical  proper¬ 
ties  of  bulk  ice  to  tire  latest  review  on  this  subject 
(Petrenko  1993).  Let  us  just  remember  that  a  small 
(compared  to  a„)  value  of  crp  results  from  lattice 
polarization  arising  from  protonic  charge  carri¬ 
ers'  motion  To  prevent  the  polarization  from  in¬ 
creasing  to  infinity,  ions  and  Bjerrum  defects 
have  to  move  in  a  consistent  manner  to  keep  their 
fluxes  equal.  An  absence  (or  reduction)  o.‘  any 
noticeable  difference  between  a„  and  a0  would 
imply  either  that  such  a  correlation  in  the  motion 
of  ions  and  Bjerrum  defects  is  absent  (or  reduced) 
or  that  partial  conductivity  of  ions  is  close  to  the 
conductivity  of  Bjerrum  defects. 

As  can  be  seen  from  Figure  13,  the  frequency 
dependence  of  as  is  much  less  pronounced  than 
that  of  ou  This  is  especially  true  for  the  ice  sur¬ 
face  after  annealing,  i.c.,  tire  smoothest  surface. 
Since  some  fraction  of  the  current  flowing  along 
the  surface  inevitably  penetrates  the  bulk,  tire 
remaining  slight  frequency  dependence  can  be 
ascribed  to  these  accompanying  bulk  currents. 
The  rougher  the  surface  is,  the  greater  is  tiro 
“bulk"  contribution  of  the  "surface"  conductivi¬ 
ty,  in  excellent  conformity  with  the  results  shown 
in  Figure  13.  Thus,  we  can  tell  th.it,  at  least  in  the 
frequency  range  of  10*1  to  3x10''  Hz,  there  is  ei¬ 
ther  no  correlation  in  tire  motion  of  ions  and  Bjer- 
runr  defects  in  the  subsurface  layer  at  all  (as  is  tire 
case  in  water)  or  that  this  c  orrelation  is  weak. 

An  alternative  explanation  for  the  weak  fre¬ 
quency  dependence  of  os  can  be  that  the  cor¬ 
responding  dielectric  relaxation  time  for  tire  sub¬ 
surface  layer  p*  might  be  several  orders  of  mag¬ 
nitude  more  than  that  of  the  bulk  ru-  Therefore, 
we  would  search  for  the  Debye  dispersion  for  tire 
liquid-like  layer  in  a  very-high-frequency  range, 
where  tlu>  standard  methods  for  measuring  o  are 
not  applicable. 

Let  us  try  to  estimate  the  specific  conductivity 
of  tire  surface  layer  o  ,s  (dimensions  of  f2"'nr  '), 
proceeding  from  lire  data  in  Table  2  and  the  re¬ 
sults  for  the  liquid-like  layer  thickness  L,  given  in 
the  previous  sections. 

Oss  =  ^  ■  (6) 

Thus,  taking  L  -  1.2xl()'4  nr  (T  =  -  IOC"  [Barer  et 
al.  llJ77])  and  os=  10  111  il  *,  we  obtain 


a^lO-1  il-'m-1. 

If  we  substitute  extreme  values  of  os  for  pure 
ice  from  I0'n  up  to  3x10"®  Sir1,  we  shall  get 

ass«(10-2-3xl01)a-1m-1.  (8) 

At  the  same  temperature,  the  bulk  low-fre¬ 
quency  conductance  Oo  falls  in  the  limits  (2.5  x 
10~®  to  6.4xlO-10)  m-1  (Petrenko  1993);  i.e.,  it  is 

from  6  to  11  orders  of  magnitude  less  than  aiS 
from  eq  8.  Tire  bulk  high-frequency  conductivity, 
a„  =  1.6xl0~5  D-1  m-1,  is  also  much  less  than  a^.  If 
we  assume  that  there  is  no  correlation  in  the  mo- 
tioi .  of  ions  and  Bjerrum  defects,  then  the  charac¬ 
teristic  dielectric  relaxation  time  of  the  surface 
layer  would  be  the  Maxwell  relaxation  time  Tms 

tms  =  (UH2-3X10-9)  s  (9) 

Oss 

which  corresponds  to  tire  frequency  range 

/=— 1 - -  (5xl07-  1.6xlo“)Hz  (10) 

2jitms 

where  to  is  tire  vacuum  dielectric  permittivity 
and  £„=  3.2. 

Unfortunately,  1  ant  not  aware  of  any  attempts 
to  measure  cjs  as  a  function  of  frequency  in  tire 
specified  range,  which  could  give  more  accurate 
values  of  ctw  in  conformity  with  eq  9.  If  in  the 
future  such  measurements  are  carried  out,  then 
this  will  also  allow  more  accurate  determination 
of  tire  liquid-like  layer  thickness  L  using  eq  6. 

The  result  of  estimating  the  charge  carrier  con¬ 
centration  in  the  liquid-like  layer  using  the  fol¬ 
lowing  formula  is  quite  interesting 

/»,  =  ?»-.  (11) 

p.Ci 

If  we  take  ass  =  3x10’  U  'nr 1  and  g  =  10  7  m2  V_l 
s'1  for  ions  and  p  ~  10-8  m2  V"1  s"1  for  Bjerrum  de¬ 
tects  (Petrenko  1993),  then  an  estimation  using  eq 
1 1  for  the  n  of  ions  and  Bjerrum  defects  gives 
3xl027  nr-3  and  5xl02K  nr"3  respectively.  (We 
used  c\  =  0.62t*  and  C)  -  0.38c  [Hubnrann  1979].) 
These  are  extremely  large  values  if  we  recall  that 
the  concentration  of  water  molecules  in  ice  is 
about  3x102k  in  3. 

The  above  estimations  can  be  much  more  accu¬ 
rate  if  we  employ  for  our  calculations  the  correla¬ 
tion  time  for  rotational  motion  of  water  mole¬ 
cules  at  the  ice  surface,  found  in  NMK  experi¬ 
ments  by  Mizuno  and  Hanafuza  (1987)  (see  Fig. 
9).  In  bulk  ice,  this  time  coincides  with  the  dielec- 


trie  relaxation  time  tp  and  is  related  to  the  con¬ 
ductivity  by  the  followin;’  formula 

tp*  =■  <I>  (aH>a  +  <M  =  =  $  ^  (12) 

\  t’l  i’t  I  <M  <\t 

where  <t>  =  3.85  k$Trtxv  rlH,  -  2.76  A  is  oxygen- 
oxygen  d’stan  .  .  "ion"  and  "13"  subscripts  corre¬ 
spond  to  ionic  and  Ujerrum  defect  conductivity 
and  i’i  and  iq  to  the  electrical  charges  of  ions  and 
Bjerrum  defects.  Substituting  Tp  =  1U10  s  (from 
Fig.  8),  T  -  263  K  and  rq  =  0.38c  into  eq  12,  we 
obtain  for  the  net  conductivity  of  the  surface  layer 

as,--  11 12-'  nr'  (13) 

i.e.,  a  value  close  to  the  upper  limit  found  in  eq  8, 
The  corresponding  concentration  of  Ujerrum  de¬ 
fects  in  the  surface  layer  according  to  eq  1 1  will  be 

»[!  -  2x  102H  m'*\  (14) 

If  there  is  no  correlation  between  fluxes  of  ions 
and  Ujerrum  defects  in  the  liquid-like  layer,  we 
should  use  eq  9.  Then 

Os, «=  0.3  U'1  nr1,  it,,  =  5x  102(’  nr-1.  ( 15) 
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riflin’  14.  Temperature  dependence  of  surface 
conductivity  of  the  (0001)  plane  for  pure  and 
doped  ice.  "O-lke"  stands  for  ice  from  C-ike  pond 
near  the  Syowa  station  in  Antarctica,  "Mendenhall" 
stands  for  ice from  Mendenhall  Glacier,  Alaska  (after 
Marnoand  Nisliitntirn  197H). 


Thus,  we  can  conclude  that  the  charge  carrier 
concentration  in  the  ice  surface  exceeds  by  many 
orders  of  magnitude  that  of  the  bulk  and  in  some 
cases  might  correspond  to  a  superionic  state,  i.e., 
a  state  in  which  all  (or  nearly  all)  molecules  are 
ionized  and  mobile. 

Caranti  and  Illingworth  (1083a)  remark  in 
their  paper  that,  in  a.c.  measurements  of  real  and 
imaginary  pans  of  the  complex  surface  conduc¬ 
tivity,  the  phase  angle  was  always  small,  show¬ 
ing  surface  currents  to  be  predominantly  resis¬ 
tive.  This  is  exactly  what  must  be  the  case  for  an 
ordinary  conductivity  (without  correlation  be¬ 
tween  charge  carrier  fluxes)  at  frequencies  much 
less  than  Tp  h 

Temperature  dependencies  of  the  suriace  con¬ 
ductivities  of  pure  and  doped  ice  are  plotted  in 
Figures  14  and  15.  The  behavior  of  this  depen¬ 
dence  is  different  at  temperatures  above  and  be¬ 
low  -6"C.  It  is  quit”  possible  that  a  sharp  rise  of  a, 
at  T  ^  — 6"C  is  associated  with  an  appearance  of 
the  liquid-like  layer  on  the  ice  surface,  which  is 
detected  by  optical  ellipsometry  and  X-ray  dif¬ 
fraction  methods  (see  previous  sections).  An 
interesting  observation  is  the  difference  in  a,  cor¬ 
responding  to  different  crystallographic  surfaces 
of  ice  (Fig.  15).  The  conductivity  of  the  basal 
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ligure  IS.  Surface  conductivity  of  ice  as  a  fane 
lion  of  crystallographic  orientation  and  ternpe-'- 
ture.  lee  single  crystals  were  taken  from  the 
Mendenhall  Glacier,  Alaska  (after  M (•■■•in  1 073) 
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Figure  1 6.  Variations  of  surface  conductance  after  evacu¬ 
ation  of  pure  VMlcr for  "pure"  natural  ice  (O-ike  pond). 
Dotted  lines  indicate  expected  variations  of  the  conduc¬ 
tance  calculated  from  measured  temperature  of  the  ice 
surface  (after  Macao  and  Nishimura  1 078). 


Figure  1 7.  Typical  thickness  dependence  of  the  mea¬ 
sured  current  intensity  obtained  for  icc  films  deposited 
onto  a  glass  substrate  of  temperature  2 20  K.  Three 
different  experimental  runs  are  shown  (after  Chrzanowski 
1988). 


(0001)  surface  was  maximum.  This  result  is  in 
agreement  with  optical  measurements  by  Elbaum 
el  al.  (1093)  bul  is  opposed  to  the  results  of  optical 
ollipsometry  (where  liquid-like  layer  thickness  is 
greater  on  the  prismatic  (1010)  surface)  and  X-ray 
diffraction  (where  the  difference  between  the 
(0001)  and  (1010)  surfaces  is  negligible). 

Maeno  (1973)  and  Maeno  and  Nishimura  (197b) 
found  that  an  intense  pumping  out  of  water 
vapors  from  above  the  ice  surface  changes  the 
surface  conductivity  considerably  (I:ig.  16).  Pump¬ 
ing  out  of  water  vapors  accelerates  evaporation 
from  the  surface  and  its  cooling.  Aside  from  tins 
trivial  effect,  a  lowering  of  water  vapor  pressure 
above  the  ice  increases  tire  surface  conductivity. 
The  phenomenon  discovered  by  Maeno  shows  the 
degree  to  which  the  conductivity  of  the  surface 
layer  is  sensitive  to  physical  conditions  at  the  ice 
interface. 

An  interesting  result  was  obtained  in  the  work 
by  Chrzanowski  (1988),  who  investigated  the  d.e 
conductivity  of  thin  ice  films  deposited  onto  a 
high  resistance  dielectric  substrate  (high  quality 
glass  or  quartz  plate).  Figure  17  demonstrates  that 
the  conductance  of  an  ice  film  grows  rapidly  with 
the  film's  thickness  until  it  thickens  to  0.3-0.5  pm. 
Then,  the  rale  of  the  conductance  growth  dimin¬ 
ishes  considerably.  This  result  implies  that  the  lay¬ 


er  close  to  the  ice /dielectric  interface  has  a  much 


larger  specific  conductivity  than  the  bulk,  and  that 
the  thickness  of  this  layer  is  about  U.3-0.5  pm.  litis 
thickness  (taking  also  into  account  the  temper¬ 
ature  of  220  K)  exceeds  considerably  values  found 
in  other  works.  I  believe  that  this  quantity  can  be 
associated  only  with  the  screening  length  hi  ice, 
which  was  considered  in  detail  by  Petrenko  and 
Ryzhkin  (1984). 

Caranti  and  I.amfri  (1987)  made  a  very  interest¬ 
ing  attempt  to  determine  charge  carrier  mobilities 
in  the  ice  surface  layer,  lire  authors  used  for  this 
purpose  measurements  of  tire  Hall  effect  and  sur¬ 
face  conductivity.  Owing  to  the  very  small  mo¬ 
bilities  of  charge  carriers  hr  ice,  tire  measurements 
of  tire  Hall  effect  are  quite  complicated.  Caranti 


and  Lamfri  used  an  a.c.  setup  with  two  distinct 


frequencies  of  tire  magnetic  field  (50  Hz)  and  the 
driving  current  (107  Hz),  measuring  tire  1  lall  effect 
at  the  total  frequency  157  Hz.  Figure  18  shows  the 
temperature  dependence  of  the  carrier  mobilities 
hr  tire  ice  surface  determined  by  Caranti  axrd  Lani- 
fri.  An  unusually  large  value  of  mobility  (3x10^* 
nr2  V"1  s'1)  is  striking.  This  is  three  orders  of  mag¬ 
nitude  greater  than  the  mobility  of  ions  in  water 
and  ice  and  four  orders  greater  than  the  mobility 
of  Bjerrunr  defects  hr  ice  (Petrenko  1993).  More¬ 
over,  such  high  mobility  can  not  be  attributed  to 
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ligutc  IS.  Hall  mobility  p  (///•'  V  1  s  1 jut  157  llzasa 
function  of  temperature  for  normal  ami  rubbed  stir 
fiu  es.  The  magnetic  field  II, ]  m/s  0.38  '1'  {after  Cnranti and 
lamfri  ISS7). 


Nonetheless,  tho  result  of  Caranti  and  Zamfri 
is  indirectly  confirmed  by  earlier  measurements 
of  Bullemer  and  Riehl  (1068).  who  tried  to  mea¬ 
sure  a  bulk  Hall  effect  in  ice  and  obtained  g  = 
1.4xll)4nrV  's  1  at  7’  -  -8‘V.  Although  liulleni- 
er  and  Riehl  used  guard  rings  at  the  Hall  elec¬ 
trodes  to  eliminate  the  effect  of  surface  currents, 
they  apparently  did  not  succeed  because  they  did 
‘observe  the  Hall  effect,  which  is  forbidden  for 
bulk  ice  ill  principle  (Sokoloff  1073,  Gosar  1074). 
Hence,  Bullemer  and  Riehl  also  determined  the 
surface  Hall  mobility.  It  is  remarkable  that  their 
result  (although  they  used  d.c.  measurements)  is 
in  a  good  agreement  with  the  Caranti  and  l.amfri 
result.  Then,  with  what  is  this  extremely  high 
mobility  associated?  Perhaps  with  electrons  in 
the  ice  surface?  We  do  not  know  the  answer. 

Pinaily,  Camp  and  Halchin  (147d)  studied  the 
electrical  properties  of  ice/ice  (grain  boundaries) 
and  ice/quart/,  interfaces  in  the  temperature 
range  from  0  to  -2l)"C  and  in  the  frequency  range 
from  d.c.  to  It)7 1 1/.  They  found  that  the  widths  of 
the  "electrical"  interfaces  (ice/ice and  ice/dielec¬ 
tric),  estimated  from  the  conductivity  mea¬ 
surements,  depend  on  the  frequency.  The  higher 


Table  3.  Summary  of  estimates  of  grain  boundary  thickness  ft  (after  Camp 
and  Halchin  1979).  The  measuring  frequency  is  indicated  in  parentheses. 


Sample 

Type  o)  measurement 

Pure  no.  2 

Prom  high  frequency 

Pure  no.  .k 

Prom  capacitance  data 

Prom  d.c.  compared  to  doped 

Pure  no.  1 

l'rom  capacitance  data 

1  roped 

Prom  capacitance  data 

Ifleel we  gram  tviwilii/y  llnekne:* 

-rc 

ft  (S  x  i(](l)  UK)  A  (order  ol  magnitude) 
ft(ll>(’)  7lX)  A 

ft  (d.c  )  24(H)  A  (onler  ol  magnitude) 

ft  ( 111*’)  KHXlA 

-pc  -14  r 

6(1(1”)  I.’  |llll  .l|lin 

ft  (1117)  S)(m  <1.1  (no 


ions  at  all.  Liven  in  the  case  of  an  activation-free 
diffusion  motion  of  ions,  their  mobility  cannot 
exceed  the  value 

g-V(T2  (16) 

Ml/ 

where  v  -  frequency  of  ion's  hops  between  mol¬ 
ecules 
c  ~  charge 

u  -  length  of  a  hop  (-  It)10  m). 

To  explain  |t  «  3x  It)  4  m2  V"1  s  ',  we  have  to 
accept  that  v  =  10ir’ s  ',  which  is  much  greater 
than  all  possible  phonon  frequencies  in  ice  (as 
well  as  in  any  other  material). 


the  frequency,  the  thinner  the  surface  layer  in 
which  the  dielectric  relaxation  time  is  small 
enough  for  water  molecules  to  be  able  to  reorient. 
That  is.  the  surface  layer  is  not  homogeneous:  the 
charge  carrier  coruentiation  is  greater  near  the 
surface  and  decreases  away  from  the  surface  into 
the  bulk  (see  Table  3,  which  is  a  summary  of 
results  of  Camp  and  Halchin). 

SURl'ACL  CHARGE  AND 
SURIACE  POTENTIAL 
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The  study  of  the  surface  potential  of  ice  and 
electric  charges  at  the  ice  surface  has  quite  a  long 
history  and  continues  to  attract  the  attention  of 


numerous  scholars.  Such  interest  in  the  electrostat¬ 
ic  properties  of  die  surface  of  ice  is  stimulated  by 
both  fundamental  problems  associated  with  it  and 
certain  practical  issues.  Fundamental  questions, 
which  tire  study  of  the  surface  potential  and  tire  ice 
surface  might  help  to  resolve,  include  tire  issue  of 
tire  structure  of  the  surface  of  ice.  For  instance,  if 
all  die  molecules  at  die  ice  surface  are  oriented 
"proton-outwards,"  as  some  dreoretical  models 
suggest  (Fletcher  1968),  diis  must  lead  to  a  positive 
charge  at  the  surface  A*  and  a  positive  potential  of 
the  ice  surface  <ps.  Bodi  A*  and  cps  are  determined  by 
die  thickness  of  the  surface  layer  d  and  the  fraction 
of  the  oriented  water  molecules  contained  in  diis 
layer.  Thus,  measurement  of  A*  and  %  in  conjunc¬ 
tion  with  other  data  (d,  for  example)  might  pro¬ 
vide  valuable  information  on  the  microscopic 
structure  and  nature  of  dro  ice  surface  layer. 

The  practical  issues  stimulating  the  study  of 
electrical  properties  of  die  ice  surface,  incorporate, 
first  of  all,  the  problem  of  atmospheric  electricity 
(see  details  in  Iiobbs  1974)  as  well  as  die  problem 
of  adhesion  and  friction  of  ice  and  snow  (Petrenko 
1994). 

Before  we  proceed  to  a  review  of  experimental 
data,  let  us  consider  basic  theoretical  relationships 
defining  the  electrostatic  properties  of  the  ice  sur¬ 
face  and  the  most  frequently  employed  experi¬ 
mental  techniques  used  to  determine  of  A*  and  cps. 

Figure  19  shows  schematically  the  distribution 
of  the  electric  charge  density  p  in  ice  near  the 
surface  (.v  -  0).  The  coordinate  scale  is  chosen  de¬ 
liberately,  First,  since  ice  is  a  conductor,  in  equilib¬ 
rium,  the  electric  field  within  the  bulk  must  be 
zero,  i.e.,  £(+  00 )  =  0.  If  we  use  the  formula  for  the 


electric  field  strength  above  a  charged  surface 


£  = 


As_ 

2eo 


(18) 


we  shall  obtain  an  apparent  result 

.  +oo 

E(+o»)  ~  -A-  p(x)dx  =  0. 
2eo  „ 


(19) 


In  other  words,  die  sum  of  all  charges  in  the  sur¬ 
face  layer  of  ice  must  be  zero. 

A  very  thin  layer  ot  die  surface  charge  per  se, 
widi  die  density  A*,,  is  composed  of  ions  adsorbed 
by  die  surface  from  the  air,  electrons  entrapped  in 
die  surface  layer,  etc.  Hus  layer  is  screened  on  die 
inner  side  by  a  much  diicker  layer  of  a  space 
charge,  composed  of  mobile  charge  carriers,  i.e., 
predominantly  Djerrum  defects  and  ions.  Using 
Poisson's  equation,  wo  obtain 

JL  =  P(_XL  (20) 

dx  Cue- 


and  as  dip/dx  =  -£,  we  can  easily  find  the  relation¬ 
ships  between  p,  £  and  <p 


<21> 


r_  is  high-frequency  dielectric  permittivity  of 
ice. 

To  estimate  a  fraction  of  water  molecules  orient¬ 
ed  along  the  electric  field,  i.e.,  widi  dipole  mo¬ 
ments  oriented  normally  to  the  surface,  we  have  to 
find  the  component  of  the  polarization  vector  p 


Pillin'  1 9.  S/huv  distribution  of  elect  no  charges  in  iee  item  the  surface. 
Seales  arearbitran/. 
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that  is  determined  by  the  reorientation  of  water 
molecules  (Pm) 


Pm  =  M  PM 
\«h2 o! 


n  H20 


(22) 


where  hh?.o  and  nL  are  water  molecule  concentra¬ 
tion  in  ice  and  the  concentration  of  water  mole¬ 
cules  oriented  along  the  field  respectively.  Sin.  e 


Pm=  Eo(£s-E~) 


E= 


then 


rX 

P  (£)d% 

.  (i 


(23) 


.'ia_  = 

"1120  M  HH20 


£»-£-  \ 
EodWllHjO/ 


p(W 


(24) 


Here,  Es  is  the  static  dielectric  permittivity  ^of 
ice  and  M  is  the  molecular  dipole  moment  (M  ~ 
8.66x10 -30  °C  m  [Hobbs  1974]). 

'lliis  molecular  "ordering"  reaches  its  maxi¬ 
mum  at  the  surface  where 


somehow  to  zero  when  x—>  °°  since  the  ice  bulk  is 
disordered.  Then  the  surface  charge  As  associated 
with  this  polarization  is 


As- 


Jo 


d—dx  =  Px(  0). 
dx 


(27) 


To  find  the  other  quantity  that  can  be  mea¬ 
sured  from  an  experiment,  namely,  the  surface 
potential  cps  (using  eq  21),  we  have  to  know  how 
the  charge  p(x)  is  distributed  in  ice;  in  other 
words,  we  have  to  know  tire  laws  of  screening  of 
an  electric  charge  in  ice.  Tire  screening  of  the 
charge  As  will  occur  differently,  depending  on  the 
ratio  of  the  average  thermal  energy  kuT  to  the 
electrostatic  energy  of  the  charge  carriers  I  c,  (ps  I . 
If  tiro  surface  potential  is  small,  so  that 


k0T  »  I  a  (ps  I  (28) 

then  in  this  case,  as  was  shown  by  Petrenko  and 
Ryzhkin  (1984) 


/14  _  ~ 

»1I20  CwMmUjO 


(^)  C(x)  =  -bi—  exp(-.xn  x)  +  exp(-tC2X)  (29) 
2eoe~  2e.ue.s 


Asa 


a 

0 


(26) 


Note  that  the  fraction  of  ordered  molecules  in 
eq  25  depends  only  upon  the  density  of  charge  Asa 
adsorbed  by  the  ice  surface.  In  other  words,  this 
result  will  still  hold  if,  owing  to  the  change  of 
temperature  or  doping,  the  concentration  of 
charge  carriers  changes,  resulting  in  a  change  of 
the  screening  charge  p(.v)  in  tire  bulk.  Therefore, 
eq  25  provides  a  remarkable  and  simple  way  to 
study  the  ordering  of  water  molecules  in  the  sur¬ 
face  layer  through  a  measurement  of  the  electric 
charge  A^,  adsorbed  in  a  thin  surface  layer.  The 
thickness  of  this  layer  is  much  less  than  the 
screening  radius  or1. 

Above  we  have  considered  the  electric  charge 
absorbed  by  (he  ice  surface  as  a  reason  for  the 
ordering  orientation  of  water  molecules.  There 
can  be  an  inverse  phenomenon  when  an  ordering 
of  water  molecules  on  the  ice  surface  generates  a 
surface  charge.  In  its  turn,  the  ordering  of  dipole 
moments  of  water  molecules  may  be  caused,  for 
example,  by  dipole-dipole  and  dipole-quadru- 
pelo  interactions  (Fletcher  1968).  Suppose  this  or¬ 
dering  is  associated  with  polarization  I’(x)  that 
reaches  its  maximum  at  the  surface  and  decays 


where  x  /  and  x'2  are  the  two  screening  lengths 
determined  by  the  major  charge  carriers  (Bjcrruin 
defects)  and  minor  charge  carriers  (ions)  respec¬ 
tively 


-1  /  eoe-JtoT 

<*-*  1  =  4/  — - 

(30) 

V  +  /14) 

*2  =  J - 

(31) 

V  (,ztl>(/ii  + 112) 

0)  ---  3.85  kuT  too. 

(32) 

Since  «  tv,,  a  large  fraction  of  electric  field 
(and  potential)  decays  according  to  eq  29  at  the 
distance  a?  j  <<  a.'  2  ■  This  happens  because  of 
screening  by  major  charge  carriers.  The  rest  of  the 
potential  decays  at  the  characteristic  length  a'2, 
determined  by  minor  charge  carriers. 

In  the  case  of  strong  fields  (large  densities  of 
the  surface  charge),  when 

k|tT  «  I  cf(ps  I  (33) 

the  screening  is  described  by  a  system  of  nonlin¬ 
ear  differential  equations,  which  apparently  does 


not  have  an  analytical  solution.  Petrenko  and 
Maeno  (1987)  have  found  a  solution  by  making 
an  approximation  that  ignores  the  difference  be¬ 
tween  fractional  values  of  the  electric  charges  of 
Bjerrum  defects  and  ions.  Assuming  C3  -  e\-  q, 
we  obtained 

£  =  a - 1 -  (34) 

2eoe«  1  +  j-g-1 

where  the  screening  length  ae-1  now  depends  on 
Xsa 

£e-,  =  4kBTeoe-  (35) 

q  Xsa 

Equation  33  can  be  used  in  the  subsurface  region 
to  a  depth  where  the  potential  I  <ps  C\  I  becomes 
comparable  to  k^T.  Further  down,  we  should  use 
eq  29. 

Before  proceeding  to  the  discussion  of  experi¬ 
mental  data,  I  shall  point  out  an  apparent 
relationship  that  exists  between  X^,  and  the  com¬ 
ponent  of  the  surface  conductivity  os;  which  io 
determined  by  the  layer  Xs  of  the  screening 
charge  of  the  protonic  carriers.  If  in  the  bulk  ( x  = 
+»o),  the  charge  carrier  concentrations  are  denot¬ 
ed  as  «/(+<*>),  then  the  corresponding  values  of 
concentrations  at  the  surface  will  be 

Hi(0)s  exp|-AJAj.  (36) 

In  reality  one  of  tire  four  types  of  charge  carri¬ 
ers  always  dominates,  hence  the  high-frequency 
surface  conductivity  is  given  by 

as  =  X  M  hi  f  »<  d*  =  Wl-b  •  (37) 

i=l  Jo 

Petrenko  (1994)  and  Petrenko  and  Colbeck  (in 
press)  studied  an  electrification  arising  from  fric¬ 
tion  between  metal  and  dielectric  sliders  on  ice. 
Thin  films  rubbing  against  the  surface  of  a  rotat¬ 
ing  ice  cylinder  were  used  as  sliders.  The  experi¬ 
ments  were  conducted  in  tire  -4.5  to-3i.5°C  tem¬ 
perature  range  and  in  the  0.5-  to  8-m  s_1  velocity 
range;  the  maximum  pressure  exerted  on  the  ice 
by  a  slider  was  1  kPa.  In  these  experiments,  we 
found  that  during  sliding,  the  sliders  pick  up  a 
significant  positive  charge  from  ice.  In  the  case  of 
metal  sliders  (stainless  steel  and  aluminum),  a 
large  potential  difference  of  1600  V  arises  be 
tween  the  ice  and  the  slider  (velocity  =  8  nr/ s,T  = 
-31.5°C),  When  dielectric  sliders  were  used  (un¬ 
der  the  same  conditions)  strong  electric  fields  £  = 


2.1  xlO6  V  nr1  were  produced  at  the  ice/slider 
interface.  From  the  meas-urements  of  the  current  I 
passing  through  an  electrometer  used  as  a  current 
meter,  we  found  that  a  coupled  ice  and  slider 
acted  as  a  charge  generator.  Thus 

/=dQ=Xjs  (38) 

dt  ds 


where  Q  =  electric  charge 

X  =  surface  charge  density  which  the  slid¬ 
er  picks  up  from  the  ice  surface 
s  =  swept  area 
f  =  time. 

Similar  results  were  obtained  in  tire  course  of  field 
experiments  on  sliding  of  alpine  skis  on  ice.  The 
value  of  the  electric  charge  density  found  is  enor¬ 
mous:  K  -  l-6xl0'2  C  m-2.  But  even  this  value 
must  be  considered  as  the  low  limit  for  X*  since, 
first,  the  slider  is  unlikely  to  pick  up  the  entire 
charge;  second,  if  the  slider  penetrated  to  some 
depth  into  the  ice,  some  fraction  of  the  unwanted 
screening  charge  can  be  entrapped  along  with  X,. 

It  is  interesting  to  make  some  estimates  based 
on  the  value  of  X*  that  we  found.  First,  such  a  high 
density  of  the  surface  charge  corresponds  to  0.01 
electron  charges  per  water  molecule  in  tire  ice  sur¬ 
face.  Second,  let  us  estimate  the  fraction  of  water 
molecules  on  the  surface  that  are  directed  per¬ 
pendicularly  to  the  surface 


'ii  no  M 


(39) 


Substituting  X*  =  1.6xl0-2  C  rrr2  into  eq  27,  we 
find  from  eq  39  that 

dii_=0.2.  (40) 

hh2o 


And  remember  this  is  the  lower  limit!  It  means 
that  a  significant  fraction  of  water  molecules  at  ice 
surface  must  be  polarized. 

Finally,  using  eq  37  to  calculate  tire  surface  con¬ 
ductivity  as,  we  obtain 

os  -  1 .6  x  10'2  C  m~2  1 .8  x  10"8  m2  V_I  s-1 
=  3xl0~1l)  £2~>  (41) 


where  the  value  of  the  mobility  of  L-defects 
(which  must  be  tire  major  charge  carriers  screen¬ 
ing  a  positively  charged  surface)  at  T  =  -10°C  (Jac- 
card  1959)  was  used.  The  estimated  value  of  crs 
found  in  eq  41  agrees  in  tire  order  of  magnitude  to 
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experimental  data  at  T  =  -10°C  (see  Fig.  14  and  15). 
Although  it  seems  unlikely  that  in  such  a  way  all 
properties  of  os  can  be  explained  (the  temperature 
dependence  of  p.4  and  as  is  different,  for  instance), 
the  conformity  found  is  hardly  purely  incidental 
and  implies  that  this  process  significantly  contrib¬ 
utes  to  the  surface  conductivity. 

Buser  and  Jaccard  (1978)  studied  a  charge  ex¬ 
change  between  very  small  (10-pm)  ice  particles 
and  various  metal  targets,  doing  the  experiments 
at  T  =  -45°C  and  collision  velocities  of  10  m/s. 
They  found  that  the  amount  and  the  sign  of  the 
charge  transferred  during  a  collision  depended 
upon  the  work  function  of  the  metal  target  as  well 
as  the  value  and  the  sign  of  the  external  electric 
field,  in  which  the  particles  were  polarized  prior  to 
collision.  To  interpret  their  experimental  results, 
Buser  and  Jaccard  assumed  the  presence  of  elec¬ 
tronic  states  at  the  ice  surface  in  an  analogy  to 
semiconductors.  These  states  participate  in  the 
electron  exchange  between  ice  and  metal  during  a 
short  (10-7  seconds,  in  tire  de scribed  experiments) 
collision.  The  lower  limit  of  the  density  of  electron 
surface  states  was  found  to  be  X^/e  £  1016  m2.  If 
we  assume  that  the  amount  of  charge  A Q  that 
passes  between  the  ice  and  the  metal  during  the 
collision  is  proportional  to  the  difference  in  their 
work  functions  (W,  -  Wm),  then  W,  can  be  defined 
as  the  point  of  intersection  with  the  horizontal  axis 
on  the  graph  of  A Q  versus  Wm.  In  this  manner, 
Buser  and  Jaccard  found  that  W;  =  4.3  eV.  A  simi¬ 
lar  result  was  obtained  in  experiments  of  the  same 
kind  by  Buser  and  Aufdcrmaur  (1977).  Caranti  et 
a).  (1985),  who  used  the  same  method  to  de¬ 
termine  the  work  function  of  ice,  found  a  some¬ 
what  smaller  value  for  W,of4.11±0.2eV. 

Another  method  used  by  many  authors  is 
Kelvin's  method.  Figure  20  illustrates  the  princi¬ 
ple  of  measurement  of  the  difference  between  a 
work  function  of  the  studied  material  (in  our  case 
14/ v  the  work  function  of  ice)  and  that  of  an  oscillat¬ 
ing  electrode  Wm.  In  this  method,  ice  (or  another 
material)  is  connected  electrically  to  the  electrode 
that  is  oscillating  above  its  surface.  If  Wm  *  W;, 
electrons  flow  from  the  ice  into  the  metal  if  Wm  > 


Vibra'ing  electrode 

\  DC  power  supply 

\  [1/ 

\Lt  % 

Lock-in 

\  lee 

amplifier 

Wj  or  in  the  opposite  direction  if  Wm  <  Wj,  until 
the  equilibrium  potential  difference  Acp  is  estab¬ 
lished 

A(f>  e  -  Wm  -  Wj .  (42) 

This  potential  difference  produces  an  electric 
field  in  the  space  between  the  ice  and  the  oscillat¬ 
ing  electrode,  which  forms  a  flat  capacitor.  When 
the  electrode  oscillates,  the  capacitance  also  os¬ 
cillates,  resulting  in  an  alternating  electric  charge 
at  the  capacitor  and  an  electric  current  in  the  cir¬ 
cuit.  At  the  voltage  V  =  -A<p,  the  magnitude  of  the 
current  is  minimal,  which  allows  Wm  -  W,  to  be 
determined.  To  find  the  absolute  value  of  Wj,  we 
need  to  find  Wm  using  some  other  method  (from 
the  extrinsic  photo-effect,  which  has  been  done 
already  for  all  pure  metals,  for  instance).  Kelvin's 
method  was  used  to  determine  the  work  function 
oficeby  Takahashi  (1969a,b,  1970),  Mazzedaetal. 
(1976),  Caranti  and  Illingworth  (1980,  1983a,b,c) 
and  Caranti  et  al.  (1985). 

TFi-  approximate  value  of  14/ j  found  in  these 
experiments  appeared  to  be  quite  sensitive  to 
such  ambient  conditions  as  humidity,  tempera  ¬ 
ture  and,  to  a  smaller  degree,  ice  doping.  A  char¬ 
acteristic  maximum  value  by  which  Wj  was  al¬ 
tered  due  to  a  change  in  the  ambient  conditions 
was  1 8Wj  I  <  200  meV.  Tims,  for  example,  Taka¬ 
hashi  (1970)  observed  a  buildup  of  a  positive  po¬ 
tential  (30-80  mV )  at  the  ice  surface  at  the  begin¬ 
ning  of  the  condensation  process.  During  further 
condensation,  the  surface  potential  changed  its 
sign  and  readied  a  value  of  -250  mV.  The  effect 
was  more  pronounced  al  die  basal  (0001)  plane  as 
compared  to  the  prismatic  plane.  Condensation 
and  evaporation  at  the  ice  surface  were  produced 
through  its  he^Fug  or  cooling,  respectively;  the 
ambient  temptmiure  was  maintained  constant. 

The  same  author  observed  a  negative  potential 
buildup  at  the  ice  surface  after  its  rubbing  (Taka¬ 
hashi  1969b).  This  result  is  in  qualitative  agree¬ 
ment  with  the  observations  of  Petrenko  and  Col- 
beck  (in  press),  although  we  found  a  much  larger 
effect. 

Mazzeda  et  al.  (1976)  studied  variance  in  Wj  as 
a  function  of  temperature.  In  these  experiments, 
the  ice  surface  was  in  equilibrium  with  the  satu¬ 
rated  vapor  throughout  the  entire  temperature 
range  of  223-273  K;  thus,  the  effect  of  condensa¬ 
tion  and  evaporation  should  not  have  manifest¬ 
ed.  Figure  21  shows  a  typical  temperature  depen¬ 
dence  of  AWim 


Figure  20.  Kelvin's  method. 


Figure  21.  Temperature  dependence  o/AW,»,  =  1  /|ej 
(W,„  -Wj)  (after  Mazzeda  et  al.  1976). 


I'igure  22.  Change  in  ice  surface  potential  when  the 
surface  was  rimed  with  supercooled  droplets  (after 
Caranti  and  Illingworth  1983c). 

AWim=I(Wm-  Wi).  (43) 

Id 

The  arrow  in  Figure  21  indicates  the  tempera¬ 
ture  (=  -30°C)  at  which  there  is  a  transition  from  a 
weak  linear  dependence  to  a  stronger  logarithmic 
dependence.  The  surface  potential  grew  more 
positive  relative  to  the  bulk  and  reference  elec¬ 
trode  when  it  approached  the  melting  point, 
which  can  be  explained  in  terms  of  the  formation 
of  a  double  charge  layer  near  die  surface,  as 
shown  in  Figure  19.  Such  charge  distribution  aris¬ 
es  if  the  molecules  at  the  surface  are  oriented 
oxygen-outwards,  entrapping  a  positive  com¬ 
pensating  charge  from  the  atmosphere. 

Caranti  and  Illingworth  (1980, 1983a, b,c)  stud¬ 
ied  variations  in  the  surface  potential  when  the 


ice  surface  was  rimed  with  supercooled  droplets 
from  an  ultrasonic  generator  (of  50-gm  radius;  1- 
to  2-g  nr3  water  content).  A  drastic  change  in  the 
surface  potential  was  observed  (Fig,  22).  The  ef¬ 
fect  of  doping  by  NH3,  NaCl  and  HF  (up  to  con¬ 
centrations  of  10-2  M)  on  the  amount  of  variation 
in  potential  was  insignificant.  Tire  authors  came 
to  the  conclusion  that  the  variation  in  the  surface 
potential  of  ice  during  its  riming  is  a  result  of  the 
disordered  ice  produced  by  rapid  freezing  of  the 
supercooled  water,  since  this  change  is  propor¬ 
tional  to  the  rate  of  freezing  and  relaxes  with  a 
time  constant  comparable  to  that  of  a  mechanical 
relaxation. 

SURFACE  ELECTROCHEMICAL 
POTENTIAL 

When  tire  impurity  concentration  in  ice  varies 
from  place  to  place,  a  static  electric  field  arises 
from  inhomogeneity  in  the  distribution  of  impu¬ 
rities,  which  are  the  effective  sources  of  protonic 
charge  carriers  in  ice.  As  a  result  of  this  inhomo¬ 
geneity  of  impurity  concentration,  a  gradient  of 
the  charge  carriers'  concentrations  n,  arises.  The 
fluxes  of  ions  and  Bjerrum  defects  are  correlated 
since  they  alter  the  orientation  of  the  water  mole¬ 
cule  when  they  pass  through  it.  The  mathemati¬ 
cal  description  of  the  general  case,  when  both  dif¬ 
fusive  and  drift  fluxes  of  all  four  types  of  charge 
carriers  must  be  taken  into  account,  is  quite  com¬ 
plicated.  Nevertheless,  we  can  easily  make  rough 
estimations  of  the  magnitude  of  the  arising  elec¬ 
tric  fields,  using  a  simplified  consideration. 

Let  us  regard  the  case  of  one  type  of  charge 
carrier  of  electric  charge  e,  and  concentration  ti,  at 
a  temperature  T.  If  a  diffusive  flux  /o  arises  from 
inhomogeneous  doping,  in  equilibrium  it  must 


be  balanced  by  the  opposite  drift  flux  /or 

,'D  =  /UK  (44) 

/ur  -  hi x  E;  Id  -  - 19,  x  grad(n,).  (45) 

Then,  employing  the  Einstein  relationship 

Dj  =  Pi  x  kyj  /  a  (46) 

we  obtain  for  the  electric  field  strength  E 

11,  x  E  1-  ( kuT/a )  x  grad(ii,).  (47) 


In  eq  44 -47,  D„  g,  and  c,  are  the  diffusion 
coefficient,  the  mobility  and  the  electric  charge  of 
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Figure  23.  Static  electric  potential  ill  the  surface  of 
freshwater  columnar  ice.  An  area  of  2. 5x2. 5  air  is 
shown.  The  columns  lie  along  the  x-axis.  The  potential 
is  given  in  volts;  T  =  -2 0°C  (from  Petrenko  1 992). 

I  he  /lit  charge  carrier  respectively;  E|t  is  Llol  Iz¬ 
mail's  constant.  For  li»i*  potential  difference  A<p 
between  two  points 

Aip  =  (cp,  -  (]J2)  =  {k,\ T/c,)  x  ln(d2l/ n ,,)  (45) 

where  ip,  and  n 2,  are  the  concentration  at  points  ! 
and  2.  Equation  48  gives  A(p  =  47  mV  for  Bjerrum 
defects  (L-and  D-defects)  at  7  =  263  K  and  »2,/h ,, 
=  It).  A  typical  magnitude  of  the  electric  field 
associated  with  Bjerrum  defects  in  ice  is  about  1- 
5  V/em  (see  Eig.  23).  Since  Aip  depends  upon  a 
ratio  of  impurity  concentration,  such  an  electric 
field  exists  even  in  very  pure  single  crystals  of  ia 
I  (Petrenko  1442)  have  studied  electric  poten¬ 
tial  on  the  surface  of  natural  and  artificial  ice 
using  a  scanning  microprobe  connected  to  an 
eleclrometi  (Eig.  24).  I  he  central  element  of  the 
experimeni.il  setup  is  an  X-Y  pen  recorder,  in 
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which  a  pen  is  replaced  by  a  thin  stainless  steel 
electrode  with  a  ball  of  0.3  mm  diameter  at  its  end. 
A  computer  controls  the  position  of  the  electrode 
scanning  over  a  given  area  of  an  ice  sample.  Scan¬ 
ning  was  from  lxl  mm  to  5  x5  cm,  with  resolution 
better  than  or  equal  to  ±0.5  mm.  The  pressure  on 
the  electrode  was  adjusted  so  finely  that  the  elec¬ 
trode  did  not  produce  visible  tracks  on  or  damage 
to  the  ice  surface.  The  velocity  of  the  mobile  elec¬ 
trode  and  tire  measuring  time  can  be  changed 
within  a  wide  range,  so  that  the  electrode  is  able  to 
measure  a  correct  value  of  an  electrochemical  sur¬ 
face  potential.  1  found  that  an  electrode  velocity  of 
less  than  10  mm  s  1  and  a  measuring  time  longer 
than  10  ms  were  satisfactory.  The  potential  was 
amplified  by  an  electrometer  preamplifier  with 
input  impedance  >  10l2ii  and  then  stored  in 
the  computer.  The  smooth  surface  of  the  ice  sam¬ 
ples  was  prepared  using  a  microtome. 

An  example  of  the  electric  potential  over  the 
surface  of  freshwater  columnar  ice  is  shown  in 
Figure  23.  The  ice  columns  lie  along  the  .v-axis.  A 
comparison  of  the  electric  potential  pattern  with 
an  optical  image  of  the  ice  sample  placed  between 
two  crossed  polarizers  reveals  their  similarity  or 
even  identity.  In  Figure  23  "rough"  and  "smooth" 
bands  belong  to  different  individual  columns.  The 
most  significant  gradient  of  potential,  and  hence 
electric  field  strength  /,'  =  grad(  ip),  is  located  at  or 
near  grain  boundaries.  Generally,  the  distribution 
of  the  electric  potential  in  ice  follows  the  structural 
pattern  of  the  ice:  arrangement  of  grains,  columns 
and  boundaries  between  the  regions  of  doped  and 
normal  ice. 

Figure  25  depicts  the  electrical  potential  over 
the  surface  of  a  two-part  ice  sample.  One  part  was 
grown  from  very  pure  deionized  water  and  one 
was  grown  from  water  doped  with  1 1) 1K  molecules 
cm  1  of  KOH.  Those  two  halves,  with  ideal  flat 
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Figure  25.  Electric  potential  on  the  surface  of  ice  sample 
assembled  of  pure  and  KOH-doped  ice.  The  steep  step  of 
the  potential  is  located  on  the  boundary  between  the  pure 
and  doped  ice.  A  "ditch"  seen  in  the  potential  of  the  doped 
(more  positive)  half  corresponds  to  locations  of  grain 
boundaries.  Dimension  of  the  area  shown  is  lxl  cm;  T  = 
-19.5’C (from  Petrenko  1992). 

and  smooth  surfaces,  were  pressed  together  for 
10  minutes  at  T  =  -5”C,  so  that  they  built  a  mono¬ 
lithic  biock.  Doping  ice  with  KOI!  provides  nega¬ 
tively  charged,  extrinsic  mobile  defects:  L-defects 
and  OH"  ions  (Howe  and  Whitworth  1989).  Their 
diffusion  into  pure  ice  makes  tire  pure  part  of  the 
sample  negatively  charged  and  the  doped  part 
positively  charged.  That  is  precisely  what  we  can 
see  in  Figure  25.  At  temperatures  higher  than 
approximately  -15°C  for  pure  ice  and  -22°C  for 
saline  ice,  the  existence  of  a  conductive,  quasi¬ 
liquid  layer  on  the  ice  surface  masks  (or  screens) 
the  electric  potential  measured  with  this  tech¬ 
nique. 

PHOTO-EMISSION  OF  ELECTRONS 
FROM  THE  ICE  SURFACE 

Nason  and  Fletcher  (1975)  studied  photo¬ 
emission  of  electrons  (extrinsic  photo-effect)  from 
the  surface  of  ice  and  water.  They  found  that  the 
threshold  of  photo-emission  for  both  ice  near 
-1G°C  and  water  near  10°C  is  about  195±5  nm. 
This  wavelength  corresponds  to  tire  energy  of  6.3 
±0.2  eV.  The  emission  efficiency  for  liquid  water 
was  about  half  that  for  ice.  Photo-emission  from 
ice  decreases  steadily  with  decreasing  tempera¬ 
ture,  apparently  without  any  shift  in  the  thresh¬ 
old  value,  and  vanishes  below  approximately 
-160°C.  Ihis  decrease  of  photo-emission  with 
temperature  was  not  related  to  the  decrease  of 
conductivity  of  ice,  since  the  same  temperature 
dependence  of  photo-emission  was  found  for  ice 
heavily  doped  with  HF  (5%).  On  the  basis  of 


similarity  of  photo-emissions  from  ice  and  water, 
the  authors  concluded  that  a  quasi-liquid  layer 
exists  at  the  ice  surface  at  temperatures  T  > 
-160°C.  Since  the  threshold  is  extremely  sensitive 
to  the  value  of  the  surface  potential,  its  weak 
sensitivity  to  temperature  variations  implies,  as 
the  authors  believe,  that  the  electric  structure  (X^, 
the  screening  length)  of  the  surface  layer  of  ice 
and  water  remains  unchanged  within  the  whole 
temperature  range  (-160  to  10°C).  Nason  and 
Fletcher  believe  that  the  most  probable  sourre  of 
photo-electrons  is  tire  process  of  ionization  of 
OH"  ions  in  the  surface  layer  of  ice  and  water. 
This  hypothesis,  obviously,  contradicts  the  fact 
that  a  change  in  pH  of  ice  from  1  to  14  does  not 
affect  photo-emission  significantly. 

At  the  present  time,  we  can  suggest  an  inter¬ 
pretation  ihat  is  different.  It  is  known  that  illumi¬ 
nating  ice  whh  light  having  a  wavelength  X  £  190 
nm  excites  a  bulk  photo-conductivity  according 
the  following  reaction 

H;0  +  lrv  ■  >  H20*  ->  H20*  +  e  (49) 

where  H2O4  denotes  an  excited  water  molecule 
(Khusnatdinov  and  Petrenko  1992).  The  elections 
produced  in  this  reaction  can  tunnel  out  of  the 
ice,  provided  that  its  surface  has  a  negative 
potential.  Since  tire  reaction  (eq  49)  begins  with 
tire  excitation  of  a  water  molecule,  it  is  not  sur¬ 
prising  that  it  is  insensitive  to  the  change  in  pH. 
The  temperature  dependence  of  photo-emission 
(as  well  as  tire  difference  in  its  values  for  ice  and 
water)  might  result  from  the  change  in  the  sur¬ 
face  electric  field,  since  this  field  causes  tunneling 
of  electrons  out  of  the  subsurface  layer  of  tire  ice. 

SURFACE  OPTICAL  ABSORPTION 
IN  INFRARED  REGIONS 

Interesting  results  about  tire  molecular  struc¬ 
ture  of  tire  ice  surface  were  obtained  most  recent¬ 
ly  from  the  study  of  infrared  optical  absorption 
spectra  (Bush  and  Devlin  1991,  Rowland  and 
Deelin  1991,  Callen  et  al.  1992,  Devlin  1992,  Nix¬ 
on  et  al.  1992,  Rowland  et  al.  1993).  Those  studies 
and  the  resulting  implications  were  supported  by 
computer  simulations  using  molecular  dynamics 
simulation  methods  (Bush  and  Devlin  1991). 

The  above-mentioned  measurements  are 
based  on  the  possibility  of  distinguishing  be¬ 
tween  lines  in  the  IR  spectra  corresponding  to  the 
stretching  and  bending  mode  oscillations  of  non- 
saturated  hydrogen  bonds  (i.e.,  dangling  bonds) 
at  the  ice  surface.  Tire  frequency  of  these  oscilla- 
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tions  appears  to  be  extremely  sensitive  to  the  way 
in  which  c  given  water  molecule  is  bonded  to 
other  molecules — by  two  or  by  three  hydrogen 
bonds  (2-coordinated  or  3  coordinated  water 
molecules).  Molecules  in  the  bulk  are  4-coordin- 
ated.  Therefore,  a  comparison  between  the  inten¬ 
sities  of  2-coordinatcd  and  3- coordinated  mole¬ 
cules  allows  us  to  estimate  their  relative  concen¬ 
trations.  The  IR  absorption  spectra  observed  in 
the  experiments  are  better  interpreted  because  of 
llie  high  sensitivity  of  the  eigen  frequencies  of 
dangling  bond  oscillations  to  tire  isotopic  substi¬ 
tution  of  K  with  U.  Moreover,  these  frequencies 
are  quite  sensitive  to  the  process  of  adsorption  of 
various  gases.  Thus,  such  measurements  allow  us 
to  study  ion  adsorption  by  the  ice  surface  at  low 
temperatures. 

Specimens  of  ice  with  large  surface-to-volume 
ratios  were  produced  either  in  the  form  of  micro¬ 
scopic  ice  clusters  suspended  in  gas  or  by  vapor 
condensation  onto  a  substrate,  forming  tlun  ice 
films.  In  experiments  by  Rowland  et  al.  (1993),  for 
example,  crystalline  ice  clusters  -25  nm  in  diame¬ 
ter  formed  when  mixtures  of  an  inert  carrier  gas 
containing  a  small  amount  of  water  vapor,  with  a 
gas/H20  ratio  of -200/1,  were  introduced  rapid¬ 
ly  into  a  cooled  cluster  cell  to  pressures  of  33  kl’a 
(V3  bar)  at  temperatures  below  -120  K.  Increas¬ 
ing  temperature  and  load  pressures  create  larger 
clusters  with  consequently  smaller  surface-to- 
volume  ratios  and  reduced  relative  surface-mode 
band  strengths  (Rowiand  et  al.  1991). 

Computer  simulations  using  molecular  dy¬ 
namics  methods  have  demonstrated  that  when 
water  molecules  collect  to  produce  a  small  clus¬ 
ter,  the  following  phenomena  occur:  1)  tire  re  are 
no  molecules  possessing  two  dangling  OH 
bonds;  2)  both  2-  and  3-coordinated  molecules 
form  on  the  surface  in  approximately  equal 
amounts,  provided  that  the  temperature  is  quite 
low  (-20  K);  and  3)  in  the  annealed  clusters,  the  3- 
coordinated  species  are  dominant.  The  concen¬ 
tration  of  dangling  bonds  is  small  («  1)  even  at 
these  low  temperatures.  The  results  of  spectral 
studies  performed  in  the  above  works  are  in 
agreement  with  the  results  of  computer  simula¬ 
tions. 

REGELATION  (REFREEZING), 

SINTERING  AND  ADHESION 

In  retrospect,  the  studies  of  the  quasi-liquid 
layer  at  the  ice  surface  began  with  studies  of  ther¬ 
mo-mechanical  effects:  regelation,  sintering  and 


adhesion.  The  data  on  the  physical  parameters  of 
the  ice  surface  layer,  obtained  during  those  stud¬ 
ies,  are  ambiguous  and  much  less  accurate  than 
those  obtained  by  modern  methods.  But  we 
should  keep  in  mind  that,  first,  it  was  these  phe¬ 
nomena  that  aroused  such  a  great  interest  in  the 
ice  surface.  Second,  these  phenomena  are  still  of 
great  practical  importance.  Therefore,  a  brief  ret¬ 
rospective  view  of  these  studies  is  presented  be¬ 
low. 

Faraday  (1850)  observed  an  interesting  phe¬ 
nomenon:  if  two  pieces  of  ice  are  brought  into 
slight  contact  at0°C  or  a  higher  temperature,  they 
freeze  together.  This  discovery  initiated  a  long 
and  interesting  argument  between  Faraday  and 
Tyndall  on  one  side  and  J.  Thomson,  his  brother 
W.  Thomson  (Lord  Kelvin)  and  Helmholtz  on  tire 
other  side.  Tire  Thomsons  and  Helmholtz  held  to 
the  opinion  that  regelation  is  caused  by  melting 
under  pressure  in  die  contact  zone  between  two 
pieces  of  ice.  When  the  pressure  is  removed,  wa¬ 
ter  freezes  again,  bonding  the  pieces  together. 
Faraday  performed  several  simple  bul  ingenious 
experiments  that  proved  the  erroneousness  of  the 
Thomsons'  hypothesis,  suggesting  a  regelation 
mechanism  that  is  very  close  to  the  modern  con¬ 
cept  of  this  process.  He  assumed  that  at  tempera¬ 
tures  close  to  lire  melting  point,  a  thin  liquid  film 
always  exists  at  the  ice  surface.  When  two  blocks 
of  ice  touch,  this  film  first  forms  a  "neck"  connect¬ 
ing  the  two  blocks.  Then,  since  the  inner  part  of 
the  neck  gradually  turns  into  a  bulk  region  of  ice, 
it  freezes,  forming  a  solid  joint  between  the  two. 

Nakaya  and  Matsuinoto  (1953)  did  some  ex¬ 
perimental  work  to  verify  the  existence  of  the 
liquid-like  layer  at  the  ice  surface.  They  suspend¬ 
ed  small  spheres  of  ice  (from  1.5  to  4  mm  in 
diameter),  prepared  from  distilled  water,  on  cot¬ 
ton  threads  in  a  pendulum-like  fashion.  One  of 
these  pendulums  was  moved  laterally,  bringing 
the  spheres  in  contact  witlr  a  very  small  force, 
estimated  by  tire  authors  to  be  less  than  0.01  dyn. 
Then,  the  suspension  point  of  one  of  the  spheres 
was  displaced  laterally  until  the  spheres  broke 
free;  the  force  can  be  calculated  from  the  value  of 
the  angle  (with  respect  to  the  vertical)  at  which 
they  separated.  Nakaya  and  Matsuinoto  also  test¬ 
ed  ice  spheres  containing  0.1%  NaCl.  In  this  case 
we  would  expect  a  thicker  film  on  the  ice  surface. 
To  explain  their  results,  the  authors  had  to  as¬ 
sume  that  a  transitional  liquid-like  layer  exists  at 
the  ice  surface  at  temperatures  below  0°C. 

Hosier  et  al.  (1957)  did  similar  but  more  accu¬ 
rately  controlled  experiments  in  a  temperature 
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Figure  26.  Separation  force  us  a  function  of  tempera¬ 
ture:  A — water-vapor  saturated  atmosphere;  B — dry 
atmosphere  (if ter  Hosier  ct  id.  1957). 


range  from  -80  to  0°C  in  botli  dr)'  and  water- 
saturated  atmospheres,  between  -80  and  -40°C, 
they  observed  no  adhesion  between  spheres  (big. 
26).  Hie  lowest  temperature  at  which  adhesion 
was  measurable  in  a  saturated  atmosphere  was 
about  -25°C,  whereas  in  a  dry  atmosphere  the  ad¬ 
hesion  practically  ceased  at  -3°C.  They  explained 
this  difference  in  terms  of  a  rapid  evaporation  of 
the  liquid-like  layer  from  the  ice  surface  in  a  dry 
air,  which  makes  it  very  thin.  Hosier  et  al.  also 
found  that  the  temperature  dependence  of  the 
force  at  which  separation  occurs  can  be  expressed 
iit  terms  of  a  Boltzmann  equation  with  Hie  activa¬ 
tion  energy  of  0.46  eV. 

Jellinek  (1957a,b;  1959;  1960;  1962),  investigat¬ 
ing  the  adhesive  properties  of  ice,  conducted 
shear  and  tensile  experiments  with  ice  frozen  on¬ 
to  stainless  steel  and  optical  flat  quartz  plates  and 
various  polymers.  In  his  tensile  experiments, 
only  cohesive  breaks  through  ice  were  observed. 
Shear  tests  resulted  in  adhesive  breaks  along  the 
interface,  the  strength  of  which  depended  on  the 
loading  rate,  surface  polishing,  temperature  and 
the  type  of  the  material  used.  A  temperature  de¬ 
pendence  of  the  shear  strength  for  the  case  of  ice- 
snow  confined  between  stainless  steel  disks  is 
shown  in  Figure  27. 

If  we  assume  that  the  shear  strength  in  shear- 
type  experiments  at  temperatures  above  ~13°C  is 
determined  by  viscous  flow  of  the  quasi-liquid 
layer,  then 


Teioc-oruturo,  °C 

Figure  27.  Strength  as  a  Junction  of  temperature  for 
snow  - ice  sandwiched  between  stainless  steel  disks,  ob¬ 
tained  by  shear.  Cross- :ectio  lull  urea  1 .54  cm2;  height  6.2 
to  0.4  cm.  Adhesive  breaks  at  T  >  -13’C  and  cohesive 
breaks  all'  <-'13  ‘C.  Fan  point  represents  an  average  of  at 
least  12  tests  (after  Jellinek  1959). 


where  q  and  L  are  film  viscosity  and  thickness 
respectively.  The  re.lio  of  q/L  found  by  Jellinek  at 
T  =  -4.5°C  was  6.99  x  104  N  s/m  for  ice-stainless 
steel  and  LSkIO4  N  s/m  for  ice-quartz,  In  eq  50 
we  have  two  unknowns,  q  and  L.  If  the  value  of  L  = 
2.7x  10-y  in,  found  by  barer  et  al.  (1977)  at  T  -  -4°C 
(see  Table  1),  is  used,  then  the  resulting  viscosity 
of  tiie  quasi-liquid  layer  between  ice  and  quartz 
will  be  q  -  0.41  Pa  s,  as  compared  to  0.017  Pa  s 
found  by  barer  et  ai.  at  oven  lower  temperatures. 

Raraty  and  Taber  (1958)  did  some  adhesion  ex¬ 
periments.  Their  results  are  essentially  in  agree¬ 
ment  with  those  of  Jellinek. 

A  large  number  of  experimental  and  theoretical 
works  dealt  with  the  issue  of  sintering  of  ice.  Sin¬ 
tering  plays  a  major  role  in  tire  process  of  snow 
densification,  which  transforms  falling  snow  into 
ice  in  mountain  glaciers  and  the  polar  caps  of 
Greenland  and  Antarctica.  Several  physical  phe¬ 
nomena,  such  as  surface  and  bulk  diffusion,  sur¬ 
face  tension,  plastic  deformation  and  water  trans¬ 
fer  through  water  vapor,  contribute  to  the  sinter¬ 
ing.  The  relative  effect  of  each  of  the  mechanisms 
above  depends  largely  on  conditions  during  sin¬ 
tering:  temperature,  humidity,  impurity  composi¬ 
tion  and  size  of  ice  particles.  Because  of  this,  inter¬ 
preting  the  data  on  ice  sintering  is  usually  very 
difficult.  The  phenomena  of  small  particle  migra- 
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lion  (Itagaki  1967a)  on  the  ice  surface  and  regela¬ 
tion  are  also  similar  to  sintering.  Studies  of  these 
tilings,  although  contributing  significantly  to  the 
development  of  practical  applications,  were  of  lit¬ 
tle  use  for  the  determination  of  the  fundamental 
parameters  of  the  ice  surface  layer  owing  to  diffi¬ 
culties  in  interpretation.  As  examples  of  excellent 
studies  on  ice  sintering,  I  can  recommend  to  inter¬ 
ested  readers  the  following,  papers:  Kingery 
(1960),  Hobbs  and  Mason  (19641),  Itagaki  (1967c) 
and  Maeno  and  libinuma  (1983). 

ADSORPTION  OT  GASES 
ON  THE  ICE  SURFACE 

When  gas  molecules  of  a  mass  in  strike  a  solid 
surface,  they  can  either  bounce  or  get  stuck  to  the 
surface  (be  entrapped)  for  a  certain  time.  The  latter 
case  is  called  an  adsorption  of  the  molecule.  If  the 
gas  molecules  are  in  thermal  equilibrium  with  the 
solid  surface,  the  equilibrium  surface  concen¬ 
tration  Cs  of  adsorbed  molecules  will  depend  on 
the  flux  of  oncoming  molecules,  i.e.,  on  the  gas 
pressure  /’  and  temperature  T  and  on  the  average 
time  that  a  molecule  spends  on  the  surface. 

Langmuir  (1916)  assumed  that  the  maximum 
concentration  of  the  adsorbed  molecules  Cm  corre¬ 
sponds  merely  to  a  mono-molecular  layer.  There¬ 
fore,  a  molecule  coming  from  a  gas  and  colliding 
with  an  already  adsorbed  molecule  must  instantly 
bounce  off  the  surface.  On  the  basis  of  this  as¬ 
sumption,  simple  kinetic  considerations  lead  to 
the  so-called  Langmuir's  adsorption  isotherm 

W  u  k  \  /’  ,.S1) 

Cm  i  h-  k i  V 

where  Vs  is  the  actual  volume  of  the  adsorbed  gas 
(which  corresponds  to  (’„)  and  I'm  (which  corre¬ 
sponds  to  Cm)  is  the  volume  ot  the  adsorbed  gas 
when  the  entire  adsorbent  surface  is  covered  with 
a  monolayer. 

A'i  =  1  x;>  .  (52) 

[2nttik\l )w'  Cm 

brunnuer  el  al.  (1938)  extended  the  result  of 
Langmuir  to  the  case  of  multilayer  adsorption. 
This  equation  is  known  as  a  UhT  equation.  They 
found  that 

I’  =  l  a(A^1)P  (53) 

V,(/'  -  /-',.)  I'm  A? 

where 


A2  =  exp 


L'(1  and  L,'t  are  the  heats  of  adsorption  (i.e.,  the 
amounts  of  heal  liberated  when  a  molecule  is  ad¬ 
sorbed)  corresponding  to  the  first  and  all  succeed¬ 
ing  layers;  l\,  is  the  saturation  pressure  of  a  gas 
over  a  plane  surface  at  temperature  T. 

As  we  can  see  from  eq  53,  a  plot  of  / ’/ l/s(P  - 1\,) 
versus  1’//’,,  must  be  a  straight  line  with  a  slope  of 
(A'2  -  1)/Vm  A2  and  an  intercept  with  a  horizontal 
axis  of  1/Vm  k\.  In  this  way  constants  Cm  and  A2 
can  be  determined;  then,  the  surface  area  of  the 
solid  can  be  calculated.  We  can  also  find  pro¬ 
vided  that  the  value  of  is  known  (it  is  usually 
taken  to  be  equal  to  the  latent  heat  01  evaporation 
of  the  adsorbate  in  a  liquid  form).  Readers  interest¬ 
ed  in  the  theories  rtf  adsorption  can  find  them  in 
the  book  by  Steele  (1974),  for  example, 

Experimental  studies  ot  gas  adsorption  on  the 
ice  surface  are  of  a  great  practical  importance  for 
the  physics  of  the  atmosphere,  meteorology,  cli¬ 
matology  and  astrophysics.  In  addition,  such  stud¬ 
ies  contribute  to  the  understanding  ol  the  ice  sur¬ 
face  structure,  since  adsorption  is  quite  sensitive  to 
both  surface  structure  and  the  presence  of  various 
defects  at  the  surface,  as  well  .is  to  the  local  fields 
(from  the  molecular  dipoles,  for  instance)  existing 
near  the  surface. 

t'.as  adsorption  by  ice  has  been  studied  in  a 
wide  temperature  range:  from  77  to  273  K,  The  ad¬ 
sorption  ol  gaseous  nitrogen  was  studied  by  Ad¬ 
amson  el  al.  (1967),  Ocampo  and  Klinger  (1983a) 
and  Schmitt  et  al.  (1987);  the  adsorption  of  C02 
was  studied  by  Ocampo  and  Klinger  (1983);  the 
adsorption  of  hydrogen  was  examined  by  Smoln- 
ehowski  (1983;  theoretical  study);  and  Orem  and 
Adamson  (1969)  studied  the  adsorption  of  non 
polar  hydrocarbon  vapors  (n-alkanes).  In  these  ex¬ 
periments,  the  surfaces  of  amorphous,  cubical  and 
hexagonal  ice.  cleaved  al  77  K  and  annealed  at 
higher  temperatures,  and  natural  snow  were  in 
vestigaled. 

C  ubic  ice,  piodueed  by  the  direct  condensation 
of  water  vapor  onto  the  walls  of  a  vessel  al  a  liquid 
nitrogen  temperature,  has  the  maximum  net  sur¬ 
face  (1 1.8  m'/'g  [Adamson  et  al.  1967]).  This  sur- 
.ace  w.is  reduced  after  the  ice  had  been  annealed, 
and  was  about  1  ncVg  in  the  case  of  natural  snow. 
1  !u  value  of  the  adsorption  heal  found  exceeds  the 
heal  of  evaporation-  27li()  128 )  mol  1  .it  77  K  for  the 
adsorption  of  N,>  reported  by  Schmitt  el  al.  (1987); 
a  semi. u  value  of  25121250  ]  /  mol  was  reported  by 
Adamson  el  al.  ( 1967}.  The  a-  Iso;  pi  ion  heat  is  gen 
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orally  reduced  after  annealing  (down  to  1067  ] 
mol-1  after  4  days  [Ocampo  and  Klinger  1983a]). 

The  results  of  studies  performed  by  Orem  and 
Adamson  (1969)  showed  that  there  exists  a  signif¬ 
icant  difference  in  the  character  of  adsorbent- 
adsorbate  interaction  above  and  below  -35°C.  Bo- 
iween  -35  and  -196°C,  ice  acts  as  a  low-energy 
adsorbent  end  the  adsorption  energy  is  associated 
with  dispersion  interactions,  with  no  effect  of  lo¬ 
cal  dipole  fields.  However,  at  temperatures  above 
-35°C,  the  ice  surface  acquires  an  adsorption 
capability  close  to  that  of  a  water  surface.  These 
observations  conform  to  the  model  of  a  quasi¬ 
liquid  film  at  the  ice  surface.  A  similar  tempera¬ 
ture  transformation  in  the  adsorption  process,  but 
at  a  lower  temperature  (-40°C),  was  also  noted  by 
Ocampo  and  Klinger  (1983a).  They  also  found  a 
considerable  decrease  in  the  entropy  of  the  ad¬ 
sorbed  gas  at  T  >  -40°C.  Such  an  entropy  decrease 
could  be  explained  in  terms  of  the  ordering  of 
adsorbate  molecu'es  in  a  strong  electric  field  of 
ordered  water  molecules  at  the  ice  surface.  As  we 
have  seen  previously,  this  also  conforms  to  the 
model  of  the  quasi-liquid  layer. 

SURFACE  ENERGY 

A  net  surface  free  energy  or  simply  surface 
energy  is  defined  as  the  work  required  to  increase 
the  interface  area  by  a  unit  area  (at  a  constant 
temperature).  We  shall  need,  for  a  theoretical 
analysis,  the  energies  of  the  interfaces  between 
solid  ice  (S'),  liquid  phase  (/,)  and  water  vapor  (V). 
We  denote  them  as  Yst Yvs  and  Yvl  for  solid/ 
liquid,  vapor/solid  and  vapor/liquid  interfaces 
respectively.  All  of  these  quantities  are  tempera¬ 
ture  dependent;  it  does  not  make  sense  to  give  a 
value  to  y  without  indicating  the  temperature.  The 
best  known  quantity  is  the  most  easily  measured, 
Yvt.,  which  can  be  determined  by  a  variety  of  ex¬ 
perimental  techniques  based  on  the  measurement 


of  the  surface  tension  of  water.  (These  techniques 
include  measurement  of  the  height  of  water  rise  in 
a  capillary  and  measurement  of  tire  wavelength  of 
surface  waves,  excited  by  an  ultrasonic  generator, 
for  example.)  At  the  melting  temperature,  Yvl  “ 
75.7  mj  nr2  and  increases  with  increasing  temper¬ 
ature.  The  two  other  energies  Ysl  and  Yvs  cannot 
be  measured  quite  so  easily,  and  to  determine 
them  a  certain  ingenuity  is  required. 

We  shall  not  describe  here  the  early  and  less 
accurate  works,  and  shall  proceed  to  an  interest¬ 
ing  technique  suggested  by  Ketcham  and  Hobbs 
(1969).  That  technique  was  also  used  later  by  I  lar¬ 
dy  (1977).  The  technique  is  based  on  the  measure¬ 
ment  of  equilibrium  angles  at  the  interfaces  be¬ 
tween  two  ice  seeds  and  water  vapor  or  water 
(Fig.  28).  As  can  be  seen  from  this  figure,  the  triple 
junction  in  a  grain-boundary  groove  is  subjected 
to  three  forces  of  surface  tension.  Therefore,  for 
equilibrium,  the  following  conditions  must  be  sat¬ 
isfied 


Ygi>  -  2 Yvs  cos(ilvs/2) 

(55) 

Ygt>  =  2Yi_s  cos(flLs/2) . 

(56) 

Since  there  are  three  unknowns  YKb-  Yls  and  yvs, 
we  need  a  third  equation  in  addition  to  eq  52  and 
53.  This  equation  can  be  tire  condition  of  equilibri¬ 
um  among  vapor,  liquid  and  ice  in  the  case  of 
incomplete  wetting  (0  <  <j>  <  re).  In  this  case 

Yvs  =  Yls  +  Yvl  cost]) .  (57) 

Using  this  technique  for  determination  of  the 
surface  energies,  Ketcham  and  Hobbs  (1969) 
found  that  at  0°C 

Yvs  =  (109  ±  3)  mj  nr7 

Yls  -  (33  ±  3)  mj  nr2  (58) 


Vapor 
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a.  Vapar/solid  interface. 
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b.  Sulid/litjuid  interface. 


Figure  28.  Schematic  of  grain-boundary  grooves  (after  Ketcham  and  1  iobba 
1969). 
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Y,;b  =  (65  ±  3)  mj  nr2. 

The  reader  can  find  the  data  on  the  values  of 
the  angle  <f>  in  the  papers  by  Knight  (1966,  1971). 
Ketcham  and  Hobbs  also  calculated  that  the 
anisotropy  of  Yvs  and  y^  with  respect  to  orienta¬ 
tion  is  probably  small  and  falls  within  the  experi¬ 
mental  error.  The  value  of  y(,(,  depends,  of  course, 
on  the  grain  boundary  type  and  must  tend  to  zero 
as  the  angle  of  misoricntation  between  the  seeds 
tends  to  zero.  I  lardy  (1977)  using  a  similar  tech¬ 
nique,  found  in  the  case  of  water/ice  interface  a 
more  exact  value 

Ym  =  (29.1  ±0.8)  mjnr2.  (59) 

Recently,  Makkonen  (in  press)  concluded  from 
results  of  measurements  of  <j>  (37°)  that  Yls  =  77  mj 
in'2  for  ice  at  -25°C. 

Finally,  Kuroda  and  Lacmann  (1982),  from  the 
number  of  broken  bonds  at  the  ice  surface  and  the 
value  of  heat  of  sublimation  of  8.5x10  2,1  J  mole¬ 
cule'1,  estimated  yvs  on  the  basal  (0001)  and  pris¬ 
matic  (1010)  faces  to  be 

Yvs  (0001)  =121  mj  m  2 

Tvs (1010)  =  128  mj  nr2-  (60) 

SUMMARY  OF 

MAIN  EXPERIMENTAL  RESULTS 

before  we  proceed  to  the  review  of  theoretical 
models  developed  to  explain  tire  structure  and 
physical  properties  of  the  ice  surface,  let  us  sum¬ 
marize  the  major  experimental  results  so  that  we 
have  a  starting  point  for  comparing  them  with 
theories.  The  reader  undoubtedly  already  antici¬ 
pates  the  difficulty  of  that  task.  The  range  of  dis¬ 
crepancies  in  the  results  of  different  authors  is  too 
large,  especially  when  different  experimental 
techniques  are  used. 

Nonetheless,  let  us  try  to  establish  tire  facts 
that  either  conform  to  the  majority  of  results  or  to 
the  most  reliable  of  them.  1  shall  also  try  to  point 
out  major  discrepancies. 

All  results  can  be  classified  into  several  issues: 

1.  The  presence  or  absence  at  the  ice  surface  of 
a  special  layer,  which  has  physical  properties  dif¬ 
ferent  from  those  of  the  bulk  ice. 

2.  The  layer  thickness  and  its  temperature  do¬ 
main  of  existence. 

3.  The  molecular  structure  of  this  layer. 

4.  A  comparison  between  i  measured  properties 


of  tire  ice  surface  layer  (diffusion,  electric  conduc¬ 
tivity,  refraction  coefficient,  etc.)  and  those  of 
water  and  ice. 

After  this  is  done,  we  shall  see  whether  it  is 
possible  to  assemble  some  consistent  picture  of 
the  ice  surface  layer. 

Existence  of  a  special  layer 

Hie  greatest  consistency  is  achieved  here.  Ex¬ 
perimental  papers  published  over  the  past  20 
years  never  address  the  question  of  the  existence 
of  the  special  layer.  All  authors  admit  to  its  exist¬ 
ence.  Let  us  join  them,  keeping  in  mind  an  im¬ 
pressive  listing  of  experimental  evidence  consid¬ 
ered  above. 

Thickness  of  the  layer  and 
the  domain  of  its  existence 

Since  in  all  cases,  a  strong  temperature  depen¬ 
dence  of  the  surface  properties  was  observed,  we 
cannot  discuss  the  question  of  the  layer  thickness 
without  referring  to  its  temperature  dependence. 

Figure  29  attempts  to  summarize  all  the  mea¬ 
surements  of  the  thickness  of  the  surface  ice  layer. 
In  those  cases  when  a  temperature  dependence  of 
some  property  of  this  layer  (e.g.,  surface 
conductivity)  rather  than  the  absolute  thickness 
was  measured,  this  temperature  dependence, 
normalized  by  the  value  at  some  temperature  To, 
is  plotted. 

At  first  glance,  the  results  presented  in  Figure 
29  are  awesome,  so  great  are  the  discrepancies. 


Figure  29.  Temperature  dependencies ;  of  quasi-liquid 
layer  thickness  taken  from  the  papers  cited  in  this 
review.  Capital  letters  at  the  curves  stand  for  first  letters 
of  the  authors'  names. 
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Tims,  the  layer  thicknesses  found  by  different 
authors  differ  by  one  to  two  orders  of  magnitude 
(at  -10°C,  L  ~  10  A  [Elbaum  et  al.  1993]  and  /.  « 
1000  A  [Golecki  and  Jaccard  1978],  for  instance), 

Tire  width  of  the  temperature  domain  of  exist¬ 
ence  also  differs  by  two  orders  of  magnitude  (AT 
=  1°C  [Elbaum  et  al.  1993]  and  AT  £  100°C  [Mi- 
zuno  and  Hanafusa  1987]).  Taking  into  account 
the  thoroughness  of  the  measurements  and  con¬ 
trol  of  tire  conditions  in  these  experiments,  I  feel 
that  it  is  most  unlikely  that  tire  described  discrep¬ 
ancies  are  attributable  to  errors  in  measurement 
or  some  other  error.  Cn  the  contrary,  I  believe 
that  all  measurements  were  correct,  but  absolute¬ 
ly  different  layers  on  the  surface  of  tire  same  ice  I|, 
were  measured,  hr  other  words,  tiro  surface  layer 
must  consist  of  several  layers,  each  possessing 
specific  physical  properties,  different  from  that  of 
tire  bulk. 

hr  fact,  optical  investigations,  based  on  tire 
measurement  of  intensity  of  the  light  reflected 
from  the  surface,  arc  sensitive  to  the  molecular 
density  of  the  layer.  Such  investigations  indicate 
tire  presence  of  an  optically  dense  film  at  the  ice 
surface  in  the  temperature  range  from  0  to 
approximately  -4  to  -5°C.  The  refractive  index  of 
litis  film  n=  1.330,  determined  by  Furukawa  etal. 
(1987a, b),  is  very  close  to  the  refractive  index  of 
water  at  0°C  n  -  1.333. 

The  thickness  of  this  thinnest,  "external"  dense 
film  depends  not  only  on  the  temperature,  but 
also  on  the  type  of  crystallographic  face,  as  well 
on  the  presence  or  an  absence  of  air  (Elbaum  et  al. 
1993).  In  approximately  tire  same  temperature 
interval  (-2.5  to  0°C),  Kouelri  et  al.  (1987)  ob¬ 
served  by  means  of  X-ray  diffraction  a  film  with 
no  long-range  order,  which  is  characteristic  of  a 
liquid  phase. 

Thus,  we  can  conclude  that  in  the  temperature 
range  from  0°C  to  "minus  several  degrees  Cel¬ 
sius,"  a  film  similar  in  structure  and  density  to 
water  exists  on  the  surface  of  pure  ice. 

A  layer  with  increased  conductivity  can  be  de¬ 
tected  at  the  ice  surface  at  much  lower  tempera¬ 
tures  (-25  to  -35°C).  Its  maximum  thickness  is 
limited  by  the  value  of  the  screening  length.  An 
idea  of  the  latter  can  bo  obtained  from  the  work 
by  Chr/.anowski  (1988):  L  =  0.5  pm  at  T  -  -53°C. 
Charge  carrier  concentration  and  conductivity  in¬ 
crease  within  tins  layer  as  they  approach  the  sur¬ 
face.  Estimates  made  in  the  section  on  the  surface 
charge  and  the  surface  potential  indicate  that  the 
concentration  of  ions  and  lljerrum  defects  near 
the  surface  can  reach  a  value  comparable  to  the 


molecular  concentration  (see  eq  14  and  15).  Such 
a  large  concentration  of  defects  itself  can  be  suffi¬ 
cient  reason  for  a  disordering  of  the  structure. 

A  layer  with  an  abnormally  large  amplitude  of 
thermal  vibrations  is  observed  at  even  lower  tem¬ 
peratures  («  -50°C  [Golecki  and  Jaccard  1978]). 
The  thickness  of  this  layer  at  T  =  -1°C  reaches  tire 
value  of  1000  A. 

Even  larger  is  the  temperature  domain  (T  2: 
-100“C)  of  a  layer  defected  by  means  of  NMR 
(Mizuno  and  Hanafusa  1987).  This  layer  can  be 
characterized  as  having  a  spin-lattice  interaction 
different  from  tlrat  of  tire  bulk,  or  as  a  layer  with 
increased  self-diffusion  coefficient  and  rate  of 
molecule  reorientation.  The  thickness  of  this  lay¬ 
er  at  tire  ice/ Aerosil  interface  is  1 1  A  at  T  =  -12°C 
(Barer  et  al.  1977). 

Molecular  structure  of  the  layer 

We  have  several  sources  of  information  about 
tire  surface  layer  structure.  Hirst,  as  was  shown  by 
Koudri  et  al.  (1987),  in  the  temperature  interval 
0  to  -2.5°C,  there  is  no  long-range  order  at  tire  ice 
surface,  wlridr  is  characteristic  of  liquids.  The  dis¬ 
ordering  probably  remains  at  lower  tempera¬ 
tures,  but  the  technique  used  does  not  have  the 
necessary  sensitivity  to  detect  Liris.  An  interesting 
fact  is  that  tire  thickness  of  the  "liquid"  layer, 
which  is  measured  by  X-ray  diffraction,  appears 
to  be  the  same  for  (0001)  and  (0110)  surfaces, 
which  contradicts cllipsoinotrical  measurements. 
Therefore,  we  cannot  rule  out  the  possibility  that 
in  these  investigations  different  layers  were  ob¬ 
served. 

Golecki  and  Jaccard  (1977,  1978),  who  used  a 
proton  channeling  technique,  observed  the  pres¬ 
ence  of  a  disordered  layer  down  to  -50°C.  This 
layer,  with  a  thickness  of  about  1000  A  at  -lnC, 
was  interpreted  as  having  much  greater  ampli¬ 
tude  of  tire  thermal  vibrations  of  oxygen  atoms 
rather  than  having  an  amorphous  structure. 

Additional  indirect  evidence  of  lire  disorder¬ 
ing  in  the  ice  surface  layer  comes,  as  we  have  seen 
in  the  section  on  surface  conductivity,  from  a 
comparison  of  NMR  data  with  ihe  data  on  the 
values  of  the  concentrations  and  the  mobilities  of 
charge  carriers  (see  eq  11-14).  We  found  that  tire 
concentration  of  defective  hydrogen  bonds  (Ujer- 
rum  defects)  equals  2xl(i2K  ur3  and  is  comparable 
to  tire  molecular  concentration  of  3xl028  nr  3.  Ob¬ 
viously,  in  this  case  it  is  impossible  to  preserve  an 
ordered  structure  interconnected  by  hydrogen 
bonds  corresponding  to  ice  1|,.  Moreover,  since 
tire  activation  energy  of  the  L-dofect  mobility 
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(0.23  eV)  is  close  to  that  of  t-1  <*  (0.29  eV,  Fig, 

9),  then  according  to  eq  11  such  high  Bjerrunr 
defect  concentration  and  consequent  inevitable 
disordering  will  persist  down  to  very  low  tem¬ 
peratures  (-60°C).  We  should  keep  in  mind,  how¬ 
ever,  that  since  the  amplitude  of  tire  NMR  signal 
decreases  with  decreasing  temperature,  the  thick¬ 
ness  of  this  layer  must  also  decrease.  Finally,  let 
us  remember  that  tire  results  of  experiments  on 
tire  surface  charge  uensity  and  gas  adsorption 
imply  that  a  fraction  of  water  molecules  in  tire 
subsurface  layer  (at  least  within  a  single  nrono- 
ntolecular  layer)  must  have  their  dipole  moments 
oriented  normally  to  the  surface.  Such  orientation 
is  also  an  ordering  of  some  kind. 


Physical  properties  of  the  surface  layer 
as  compared  with  water  and  ice 

As  we  discussed  earlier,  certain  physical  prop¬ 
erties  of  the  surface  ice  layer  are  intermediate 
between  water  and  ice,  This  holds  for  viscosity 
and  self-diffusion  coefficient,  for  instance.  But 
tli ere  are  two  remarkable  exceptions  to  this  rule. 
First,  the  optical  density  of  the  quasi-liquid  film 
practically  coincides  with  that  of  water  at  0°C. 
Second,  the  net  electric  conductivity  of  tire  ice 
surface  exceeds  by  many  orders  of  magnitude  the 
electric  conductivities  of  both  ice  and  water  and 
corresponds  to  tire  upper  boundary  of  conductiv¬ 
ity  of  superionic  materials. 
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PART  II.  THEORETICAL  MODELS  OF  THE  ICE  SURFACE 


We  have  seen  in  Part  I  how  complicated  the 
structure  of  the  ice  surface  is  and  what  diverse 
physical  properties  it  possesses.  It  is  clear,  there¬ 
fore,  how  difficult  it  is  to  create  a  theoretical  model 
describing  such  physical  objects  as  the  surface  lay¬ 
er  of  ice,  and  nobody  has  yet  been  able  to  devise 
one.  Theoretical  models  do  exist,  although  there 
are  very  few,  that  in  tire  best  case  can  describe 
satisfactorily  only  one  of  the  phenomena  associat¬ 
ed  with  tire  ice  surface  in  a  limited  temperature 
range.  1  lowever,  probably  all  tire  physical  reasons 
that  make  the  ice  surface  as  it  is  were  analyzed  in 
these  works.  Therefore,  it  seems  reasonable,  in 
spite  of  a  very  limited  application  of  existing  theo¬ 
ries,  to  consider  and  discuss  them  briefly. 

lire  most  simple  and  clear  is  the  phenomeno¬ 
logical  model  based  on  tire  idea  of  a  possible  gain 
in  tire  free  energy  of  tire  ice/ vapor  interface  result¬ 
ing  from  the  formation  of  a  liquid  film  on  tire  ice 
surface.  Obviously,  this  is  possible  only  if 

Tvs  >  Yi.s  +  Yvl  .  (61) 

Within  the  framework  of  this  theory,  the  sur¬ 
face  layer  of  ice  is  considered  to  be  simple,  ordi¬ 
nary  water.  If  we  denote  dre  change  in  tire  Gibbs 
free  energy  per  unit  area  resulting  front  melting  as 
AG,  obviously  in  thermodynamical  equilibrium 

-■£?-)  =  pi.Ajtfd)  -  0  (62) 

Od  ft 

where  p  is  the  chemical  potential,  p[  is  the  water 
density  and  d  is  tire  water  film  thickness.  The  wa¬ 
ter  film  thickness  d  enters  into  AG  twice 

AG  =  (pL  SL  AT)  d  )  A y(d)  (631 

where  AT  -  7’m  -  T,  S[,  is  the  entropy  of  fusion  per 
unit  mass  and  Ay  is  tire  change  in  tire  free  energy  of 
the  surfaces 


Ay{d)  =  Yvi+Ysl-YvS'  (64) 

Ay  depends  on  tire  film's  thickness  implicitly  be¬ 
cause,  at  thickness  values  comparable  to  the  inter¬ 
atomic  distances,  two  interfaces  begin  to  "sense" 
the  presence  of  each  other.  Usually,  it  is  assumed 
that 


Ay (d)  =  A  y(°°) 


(65) 


where  r0  is  a  constant  on  the  order  of  interatomic 
distance  (see  for  example  Dash  [1991]  and  Takagi 
[1990]).  A  solution  of  tire  system  of  eq  62-65  gives 
for  tire  water  film  thickness 


2rl  A  t~)  VA 
,pL  Sl  (T,n  —  T)j 


(66) 


lire  reader  can  find  examples  of  these  calcula¬ 
tions  in  tire  paper  by  Dash  (1991),  who  also  consid¬ 
ered  an  ice/solid  interface  and  took  into  account  a 
small  correction  due  to  tlrermomolecular  pressure 
acting  oir  the  melted  layer;  also  in  tire  paper  by 
Takagi  (1990),  who  incorporated  corrections  due 
to  surface  curvature.  Figure  30  shows  the  results  of 
such  a  calculation. 

Fukuta  (1987)  and  Makkonen  (in  press)  sug¬ 
gested  that  tire  surface  can  melt  because  of  the 
pressure  applied  to  the  subsurface  layer  by  tire 
bulk,  lire  nature  of  this  pressure  becomes  clear  if 
we  recall  that  a  molecule  on  the  surface  is  subject¬ 
ed  to  tire  attraction  forces  from  the  rest  of  the 
molecules  and  tire  net  force  is  directed  inside  the 
crystal.  Since  the  pressure  lowers  tire  melting  tem¬ 
perature  of  ice,  the  surface  melts.  Fukuta,  analyz¬ 
ing  the  intermolecular  forces  and  using  the  depen¬ 
dence  Tm(P),  valid  for  microscopic  quantities  of 
water  and  ice,  estimated  the  thickness  of  tire  liquid 
film  as 

d  -  -l-5jrm  .  (6 7) 

VT.7TT 


Figure  30.  Calculated  temperature  dependence  of  a  thick¬ 
ness  of  liquid  film  on  an  icc  surface  (after  Takagi  1990). 
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Makkonen  (in  press),  who  found  experimentally 
that  eq  61  is  not  valid,  predicted 

d  =  !£ .  (68) 

(Tm  -  T) 

Within  the  context  of  Fukuta's  model,  the  liquid 
film  must  vanish  at  temperatures  T  <  ~22°C, 
where  tire  pressure  just  transforms  ice  I(,  into  its 
high  pressure  forms  (in,  Im  and  Itx)  rather  than 
into  liquid, 

Thus,  the  models  just  described  could  have 
been  applied  to  the  optical  observations  of  a  dense 
liquid  film  or  the  ice  surface  close  to  the  melting 
point,  although  the  estimation  usiirg  eq  o7  gives  a 
much  thinner  thickness  compared  to  experimental 
data  (for  example  at  T  =  -1°C;  see  Fig.  29).  In  addi¬ 
tion,  these  phenomenological  models  can  give 
nothing  more  in  terms  of  propeities  of  the  liquid 
film  than  tire  properties  of  ordinary  water.  And 
these  properties  do  not  coincide  at  all  with  the 
properties  of  the  liquid-like  layer  of  ice,  as  we 
have  seen  previously  (electrical  conductivity  and 
self-diffusion  coefficient,  for  instance). 

During  a  rather  long  period,  a  model  proposed 
by  Fletcher  has  attracted  (and  still  continues  to 
attract)  considerable  attention  (see  Fletcher  1962, 
1963, 1968).  (Note  that  Fletcher  recently  returned 
to  the  problem  of  tire  ice  surface  structure  at  much 
lower  temperatures  [see  Fletcher  1992].)  The  key 
idea  in  Fletcher's  (1968)  model  is  the  assumption 
that,  at  the  ice  surface,  a  significant  fraction  of 
molecules  (denoted  as  Uo)  is  oriented  with  their 
dipole  moments  (i.e.,  protons)  pointing  outwards, 
as  opposed  to  the  bulk  where  molecules  are  ori¬ 
ented  randomly.  This  would  result  in  a  largo 
buildup  of  positive  polarization  charge  at  the  sur¬ 
face.  Fletcher  came  to  his  conclusion  about  mole¬ 
cule  orientational  ordering  by  analyzing  tire  inter¬ 
action  between  dipole  and  quadrupole  moments 
of  water  molecules  dose  to  the  surface.  Fletcher 
estimated  the  energy  of  this  interaction  as  10  20  J 
per  molecule,  which  exceeds  an  average  thermal 
energy  ky T  «  4x10  21  ]. 

As  we  have  seen  in  the  section  on  surface 
charge  and  surface  potential,  a  positive  electric 
charge  of  a  high  density,  which  is  unambiguously 
related  to  the  water  molecule  dipole  ordering,  was 
indeed  discovered  on  the  ice  surface. 

All  other  predictions  of  Fletcher's  model  are 
consequences  of  the  molecule  ordering  in  the  sur¬ 
face  layer.  In  fact,  since  in  the  bulk  there  is  no 
ordering  (a„  =  V2),  between  the  surface  and  the 
bulk,  there  should  be  a  transitional  layer  with  a 
large  concentration  of  defective  hydrogen  bonds 


allowing  molecules  to  reorient.  This  transitional 
disordered  layer,  according  to  Fietcher,  is  the  liq¬ 
uid-like  layer.  Fletcher  finds  tine  equilibrium  thick¬ 
ness  of  this  layer  thermodynamically  by  minimiz¬ 
ing  the  Gibbs  free  energy  of  the  surface  and  incor¬ 
porating  the  following  components  of  this  energy: 
the  energy  associated  with  the  surface  orientation 
of  the  layer;  the  energy  associated  with  the  interac¬ 
tion  of  tine  molecular  quadruples  with  the  gradi¬ 
ent  of  the  electric  field;  the  interaction  of  molecular 
dipoles  with  the  total  electric  field;  and  the  interfa¬ 
cial  energy  between  tine  liquid-like  layer  and  the 
ice.  In  addition,  a  "volume"  reduction  in  G  be¬ 
cause  of  the  melting  of  the  layer  and  a  lowering  of 
the  entropy  because  of  molecule  orientation  are 
apparently  also  considered.  An  unknown  fitting 
parameter  of  tine  theory  is  e2,  the  energy  required 
to  form  an  L-defeet  in  the  liquid-like  layer.  Finally, 
Fletcher  minimizes  the  free  energy  G  as  a  function 
of  a„  and  the  layer  thickness  d.  I  deliberately  do 
not  present  tine  mathematical  details  of  Fletcher's 
calculation  because  it  is  quite  cumbersome  and, 
more  importantly,  might  be  made  more  accurate 
in  the  future.  Nevertheless,  Fletcher's  analysis  of 
tine  problem  and  tine  insight  that  he  has  achieved 
into  the  physical  nature  of  the  layer  are  absolutely 
remarkable.  Let  us  note  also  that  within  the  con¬ 
text  of  this  model,  an  unusually  high  electrical 
conductivity  of  line  surface  layer  has  an  obvious 
explanation:  it  is  ascribablc  to  a  high  concentration 
of  L-defects  in  the  liquid-like  layer.  Figure  31 


Figure  31.  Calculated  dependence  of  the  equilibrium 
thickness,  L,  of  the  liquid-like  layer  on  ice  for  various 
assumed  values  ofc2,  the  energy  of  formation  of  diffuse 
L-defects  in  the  liquid-like  layer  (after  Fletcher  1968). 
a—t2  =  6.25x10  W;  'b-~e2  =  <>-0188  cV; 
c— e2  =  0.0625  eV;  d— 1:2  =  0.1875  cV. 
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Figure  32.  Schematic  representation  of  protonic  point  defects 
concentration,  n„  as  a  function  of  a  distance  from  the  ice 
surface,  x;  d  is  the  thickness  of  a  superionic  layer.  The 
axis  is  logarithmic. 


shows  the  calculated  thickness  of  the  liquid-like 
layer  as  a  function  of  temperature.  Later,  in  his 
invited  review  paper,  Fletcher  (1973)  gave  an  ap¬ 
proximate  expression  for  the  thickness  d  of  the 
liquid-like  layer 

d{k)  =  (20  to  50)  -  25  log10(273  -  T).  (69) 

Fletcher's  (1992)  recent  calculations  suggest 
that,  at  much  lower  temperatures,  the  ice  surface  is 
reconstructed  by  an  orientation  of  surface  dipoles 
into  a  superlattice.  Such  a  reconstruction  should 
occur  on  the  (0001)  basal  plane  below  approxi¬ 
mately  30  K  and  on  the  (0110)  prismatic  plane 
below  approximately  70  K. 

On  the  basis  of  the  model  of  the  superionic  state 
of  ice  invented  by  Ryzhkin  (1985),  Ossipyan  and 
Petrenko  (1988)  suggested  a  purely  electrostatic 
reason  for  surface  melting  that  accounts  for  a  high 
surface  electric  conductivity  and  the  appearance  of 
a  rather  thick  (up  to  the  screening  length)  disor¬ 
dered  layer.  The  reasons  for  the  formation  of  tire 
layer  arc  as  follows.  As  v/as  shown  by  Ryzhkin 
(1985),  upon  tire  increase  of  charge  carrier  concen¬ 
tration  above  a  certain  critical  value  (on  the  order 
of  lO41  to  10-4  of  tire  molecular  concentration),  tire 
activation  energy  of  ion  and  Bjerrum  defect  forma¬ 
tion  drastically  decreases,  leading  to  a  catastrophic 
growth  of  defect  concentration  up  to  0.1-1  of  the 
molecular  concentration.  This  is  equivalent  to  a 
phase  transition  into  a  superionic  slate,  which  can 
not  have  a  regular  lattice  owing  to  a  high  defect 
concentration. 

If  the  surface  of  ice  is  charged  (which  is  the  case 
in  reality,  as  we  have  already  seen  in  the  section  on 
surface  charge  and  surface  potential),  this  charge 
will  be  screened  by  ions  and  Bjerrum  defects  com¬ 
ing  from  the  bulk,  as  shown  in  Figure  32.  By  taking 
into  account  the  transition  into  a  disordered 


superionic  state,  we  can  transform  a  smooth  curve 
of  the  screening  charge  concentration  distribution 
into  the  curve  with  a  thin  superionic  layer  near  the 
surface.  Within  the  framework  of  this  model,  the 
high  surface  conductivity  of  ice,  the  formation  of 
the  quasi-liquid  disordered  layer  and  the  sharp 
interphase  boundary  (detectable  by  optical  meth¬ 
ods)  can  be  quantitatively  accounted  for.  To  my 
knowledge.  Dr.  Ryzhkin  is  currently  working  on 
the  improvement  of  this  model. 

There  exists  a  very  nontrivial  and  unusual  mod¬ 
el  of  the  structure  and  properties  of  the  ice  layer  at 
temperatures  close  to  the  melting  point.  Dzy- 
loshinskii  et  al.  (1961)  showed  that  a  fluid  that  is 
less  polarizable  than  its  solid  phase  is  attracted  to 
the  solid,  causing  a  liquid  film  to  grow  at  tire  sol¬ 
id/vapor  interface.  Elbaum  and  Schick  (1991)  ap¬ 
plied  th’s  idea  to  the  ice  surface.  Figure  33  shows 
the  result  of  their  calculations.  If  we  compare  this 


Figure  33.  Calculated  temperature  dependence  of  the 
liquid  film  on  ice  thickness  (after  Elbaum  and  Schick 
1991). 
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figure  to  the  experimental  results  (Fig.  29),  it  be  ¬ 
comes  apparent  fhat  the  contribution  of  this 
mechanism  to  the  experimentally  measured 
thickness  of  the  quasi-liquid  film  is  extremely 
small.  Kroes  (1992)  used  ?  molecular  dynamics 


computer  simulation  to  study  the  surface  of  a 
small  cluster  of  water  molecules.  He  found  such 
ice  surface  characteristics  as  premelting  and  par¬ 
tial  polarization  of  molecular  dipoles. 
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Time-stepping  scheme; 

Scmi-implicit  finite  difference 

Implicit  discretization  of  the  gradient  of  the  surface  elevation  in  the 
momentum  equations  and  the  velocity  in  the  vertically  integrated  continuity 
equation 

Explicit  discretization  of  the  convective,  turbulence,  and  certain  body  force 
terms 

Numerical  solution  technique:  Finite  difference 

Type  of  grid:  Cartesian  finite  difference,  constant  grid  spacing 

Coordinate  transformations: 

Horizontal:  none 
Vertical:  none 

Example  applications:  San  Francisco  Bay;  Lagoon  of  Venice 
Typical  At  for  3-D  harocl-uic  circulation:  1  to  2  minutes 
References: 

Casulli  and  Cheng  (1992) 

Cheng,  Casulli,  and  Gartner  (1993) 
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Sheng-Johnson  Estuarine  Hydrodynamic  Model 


Author(s):  Y.  Peter  Sherig  and  Billy  H.  Johnson 
Organization:  University  of  Florida,  Gainesville,  and  WES 

Contributors  (modifications  and  enhancements): 

Name  o','  model:  CH3D  (Curvilinear-Grid  Hydrodynamic  3-D)  /  CH3D-WES 
(WES  version  of  CH3D) 

Related  models: 

EHSM3D  (Estuary  Hydrodynamic  and  Sediment  Model  3-D,  Cartesian 
coordinates) 

CH2D  (Curvilinear-Grid  Hydrodynamic  2-D  Vertically  Averaged) 

Remarks:  Dr.  Shcng  continues  to  update  and  make  enhancements  to  the 
model.  He  has  developed  water  quality/scdimcnt  transport  models  that 
interface  with  CH3D  and  EHSM3D. 

WES  has  made  significant  modifications  to  the  code  that  have  resulted  in  a 
functionally  sepaiate  version,  CH3D-WES.  The  WES  version  has  been 
modified  to  use  either  a  otransform  or  a  •/- transform,  dependent  upon  applica¬ 
tion.  Welling  and  drying  is  not  implemented  in  the  current  WES  version  of 
the  code.  A  water  quality /aquatic  ecology  model  has  been  developed  to 
interface  with  CII3D-WES.  A  3-D  sediment  transport  model  is  currently 
under  development  for  WES  (anticipated  completion  date:  October  11/93). 

Dimensionality:  3-D 

Governing  equations: 

Fully  nonlinear  Reynolds  form  of  Nnvicr-Stokcs  with  hydrostatic 
approximation 

Continuity  of  Mass 

Advcclion-Dil fusion  of  Salinity 

Advection-Diflusion  of  Temperature 

Equation  of  Slate 

Turbulence  submodel: 

Horizontal:  Eddy  viscosity 
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Vertical :  Simplified  second  moment  closure  model 
Time-stepping  scheme: 

Split  mode  solution; 

Scrni-implicil  finite  difference  for  the  external  mode  (vertically  averaged) 
solution 

Ail  terms  of  the  internal  mode  solution  are  treated  explicitly,  except  for  the 
vertical  turbulent  diffusion  terms  and  the  bottom  friction  and  surface  slope 
terms  in  the  momentum  equations 

Numerical  solution  technique:  Finite  difference 

Type  of  grid:  Curvilinear,  staggered,  boundary  fitted,  structured 

Coordinate  transformations: 

Vertical;  ocoordiuate  system 

Horizontal:  curvilinear,  boundary  lilted 

Example  applications: 

Lake  Okeechobee;  Sarasota  Hay  (Sheng) 

Chesapeake  Hay:  Delaware  Hay;  Green  Uay  (WHS) 

Typical  At  for  3-D  hnrocliuic  circulation:  2  to  10  minutes  (Hqual  time-steps 
arc  used  for  internal  and  external  modes.) 

References: 

Sheng  (1993) 

Johnson  el  al.  (1991) 

Sheng,  Parker,  and  ileiin  (1986) 
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Edinger-Buchak  Estuarine  Hydrodynamic  Model 

Aulhor(s):  John  E.  Edingcr  and  Edward  M.  Budiak 

Organization:  J.  E.  Edingcr  Associates,  Inc.,  Wayne,  PA 

Contributors  (modifications  and  enhancements):  Zcnyu  Zhang,  University 
of  Pennsylvania;  Lixing  Wu,  University  of  Pennsylvania 

Name  of  model:  GLLVHT5  (Generalized,  Longiuidinal-Latcral-Vcrtical 
Hydrodynamic  and  Transport  Model,  Version  5) 

Related  models: 

GLVHT  (Generalized,  Longitudinal- Vertical  Hydrodynamics  and 
Transport),  the  hydrodynamic  basis  of  CE-QUAL-W2,  the  Corps’  laterally 
averaged  reservoir  and  estuary  model 

GLUT  (Generalized  Longitudinal  Hydrodynamics  and  Transport),  a  1-D 
river  and  estuary  model 

Remni'ks:  Edingcr  Associates’  models  are  designed  for  quick  setup  and 
moderate  detail.  According  to  Edingcr  Associates,  this  method  permits 
(a)  rapid  issue  identification;  (b)  low  computing  costs  and  moderate  hardware 
requirements  (an  ordinary  personal  computer  or  workstation  is  sufficient  for 
most  studies);  and  (c)  iterative  refinement  of  the  modeling  procedure,  All  of 
Edinger’s  models  include  concurrent  water  quality  subroutines.  Mass  balances 
within  the  defined  computational  cells  are  perfect  to  machine  accuracy. 

Wetting  and  drying  of  the  computational  grid  is  supported. 

Dimensionality:  3-D 

Governing  equations: 

Lully  nonlinear  Reynolds  form  of  the  Navicr-Slokcs  equations  with 
hydrostatic  approximation 

Continuity  of  Mass 

Constituent  Transport 

Equation  of  Stale  relating  water  density  to  salinity  and  temperature 

Turbulence  submodel: 

Horizontal:  User  may  specify  a  variety  of  forms  for  computing  Reynolds 
stresses. 
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Vertical:  User  may  specify  a  variety  of  forms  for  computing  Reynolds 
stresses,  including  Mcllor-Yamada 

Time-stepping  scheme:  A  fully  implicit  integration  e:  vides  the  time-step  into 
two  halves.  On  the  first  half-time-step,  the  x-  and  y-momenlum  equations  are 
substituted  into  the  continuity  equation  and  the  resulting  linear  equations  are 
solved  for  the  water-surface  elevation.  The  horizontal  velocities  are  then 
solved  algebraically  from  the  momentum  equations,  followed  by  the  vertical 
velocities  from  continuity.  A  second  pass  through  this  procedure  assures  a 
perfect  water  balance.  The  constituent  and  density  equations  arc  then  solved  to 
complete  the  first  half-time-step. 

On  the  second  half-time-step,  the  horizontal  velocities  are  solved  vertically 
implicitly,  then  the  vertical  velocities  and  the  water-surface  elevation  are 
solved  explicitly.  The  constituent  and  density  equations  are  solved  to  complete 
the  second  half-time-step.  If  there  arc  significant  water-surface  elevation 
changes,  a  laycr/cell  add/suhiraet  module  is  activated. 

For  very  large  grids,  the  code  has  an  option  to  compute  the  main  variables 
with  a  similar  second-order  ADI  method.  For  this  method  the  first  half-time- 
step  solves  for  a  vertically  implicitly  and  v-z  by  substituting  the  x-momentum 
equation  into  continuity,  The  second  half-time-step  reverses  this  procedure 
and  solves  foi  v  vertically  implicitly  and  u-z  simultaneously.  The  solutions  are 
obtained  using  the  tridiagonal  algorithm.  The  ADI  scheme  is  suited  to  very 
large  giids  in  which  its  inherent  biasing  of  the  buoyancy  terms  is  negligible. 

A  third  option  is  to  use  a  firsl-ordei-in-lime,  second-urder-in-space  single- 
stage  compulation.  The  lallet  option  is  used  for  fast,  preliminary  simulations, 

Numerical  solution  technique:  Finite  difference 

Type  of  grid:  Orthogonal  finite  diflcrcncc 

Coordinate  transformations: 

Vertical:  variable  Ac,  no  transformation 

Horizontal:  A.\  *  Ay,  no  transformation 

Example  applications:  Sterling  Harbor,  Texas;  Grand  Lake,  New  Brunswick; 
Patuxent  River  Estuary,  Maryland;  Bellcdunc  Harbor,  New  Brunswick; 

Neehako  Reservoir,  British  Columbia;  East  Waterway,  Even  u,  Washington; 
Farley  Storage  Pond,  Dothan,  Alabama;  Webber  Cove- Barbara  Weil  Estuary, 
Prince  Edward  Island;  Moller  Bay,  Alaska;  Delaware  Estuary,  Wilmington, 
Delaware;  Candarli  Korlc/i  (on  the  Aegean  Sea),  Turkey;  Clinton  Lake, 
Illinois;  Du  Font  Cooper  River  Plant,  Cooper  River,  South  Carolina;  Sheldon 
Point,  Saint  John  Harbour,  New  Brunswick;  Comanche  Peak,  Texas;  Barbers 
Point  Oulfall/Barbers  Point  to  Ewa  Beach  coastal  area,  Hawaii;  Maggie 
Crcok/I  lumboldl  River,  Nevada 
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Typical  At  lor  3-D  baroclinic  circulation:  7.5  to  15  minutes 
References: 

Buchak,  Edward  M.  (1993).  Personal  communication,  14  October. 
Edinger  and  Buchak  (1980,  1985) 
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Zhao-Tabios-Shen  Estuarine  Hydrodynamic 
Model 


Authors):  Di  Hua  Zhao,  Guillermo  Q.  Tabios,  Hsich  Wen  Shcn 
Organization:  University  of  California,  Berkeley 
Contributors  (modifications  and  enhancements): 

Name  of  model:  RBFVM-2D 
Related  models:  none 

Remarks:  This  modci  is  included  for  review  because  it  is  the  only  published 
finite  volume  model  currently  being  used  for  surface  water  flows.  In  its 
present  form  it  is  not  applicable  to  stratified  estuarine  How;  however,  there  is 
great  potential  for  3-D  model  development.  The  finite  volume  method  has 
merit  for  shallow  surface  water  flow  problems  because  wetting  and  drying  of 
subdomains  is  handled  naturally.  An  implicit  version  is  being  developed  to 
alleviate  the  lime-step  restriction. 

Dimensionality:  2-D  vertically  averaged  hydrodynamics 

Governing  equations:  2-D  shallow-water  equations  with  constant  density 

Turbulence  submodel: 

Horizontal:  neglected 
Vertical:  none 

Time-stepping  scheme:  Explicit 

Numerical  solution  technique:  Finite  volume  method 

Type  of  grid:  Unstructured  with  quadrilateral  and  triangle  elements 

Coordinate  transformations: 

Horizontal:  none 
Vertical:  none 

Example  applications:  Kissimmee  River  Basin 

Typical  At  for  2-D  vertically  averaged  circulation:  less  than  1  second 
Reference:  Zhao,  Tabios,  and  Shcn  (1992) 
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7  Recommended  Criteria  for 
Selecting  a  Hydrodynamic/ 
Salinity  Model  for  the 
LSJRE 


The  Slate  of  Florida  has  tasked  ihe  SJRWMD  to  determine  ihe  ecological 
and  water  quality  condition  of  the  St.  Johns  River  and  to  develop  management 
strategies  that  balance  multiple  beneficial  uses  of  the  river.  Potential  beneficial 
uses  may  include  water  supply,  navigation,  industrial  and  residential  use  of 
shorelines,  domestic  and  industrial  waste  discharges,  fisheries,  wildlife  habitat, 
and  recreation.  From  an  engineering  perspective,  the  objectives  of  the 
SJRWMD  are  to: 

a.  Understand  the  existing  physical,  chemical,  and  biological  functions  of 
the  LSJRE. 

b.  Develop  tools  to  scientifically  predict  the  impacts  of  alternative 
management  plans  and  engineering  projects  on  the  LSJRE. 

c.  Develop  a  methodology  for  evaluating  the  costs  and  benefits  of 
alternative  management  plans,  based  on  predictions. 

Numerical  models  have  ihe  potential  for  integrating  basinwide  hydrodynamic 
and  water  quality  data  with  scientific  principles  in  a  quantitative  manner. 

Given  the  directive  of  the  State  of  Florida  to  the  SJRWMD,  such  sophisticated 
analysis  tools  arc  warranted. 

The  selection  of  a  hydrodynamic  model  lor  the  LSJRE  is  the  first  step  in 
creating  a  comprehensive  modeling  and  analysis  system.  The  recommended 
features  of  a  model  that  will  meet  the  long-range  goals  of  the  SJRWMD,  and  a 
brief  justification  of  each  feature,  are  given  below: 

a.  Minimum  dimensionality:  2-D,  vertically  averaged  estuarine 
hydrodynamics. 

b.  Potential  dimensionality:  3-D,  baroclinic  hydrodynamics. 


48 


Chapter  7  Recommended  C  'teria  tor  Selecting  a  Hydrodynamic/Salinity 


t\  Geometry:  boundary-fitted  grid  with  local  refinement  options. 

d.  Welting  and  dr,’ing  option:  required. 

e.  Horizontal  turbulent  fluctuation  terms:  required 

f.  Vertical  turbulent  closure  model:  documentation  required. 

g.  Mode I  maintenance:  required 

A  brief  justification  of  each  feature  is  given  in  the  following  sections. 


Minimum  Dimensionality 

The  LSJRE  includes  approximately  110  miles  (177  km)  of  tidally  influ¬ 
enced  river  channel.  Although  several  1-D  models  have  been  applied  to  ihe 
LSJRE,  the  process  of  verifying  a  1-D  model  frequently  results  in  unrealistic 
values  for  roughness,  channel  geometry,  and  water  quality  diffusion  coeffi¬ 
cients.  Two-dimensional,  vertically  averaged  models  have  been  found  to 
require  far  less  "tuning"  to  achieve  satisfactory  tidal  simulations,  and  arc  more 
useful  for  making  valid  predictions. 

The  lower  100  miles  (160  km)  of  the  St.  Johns  River  varies  from  600  ft 
(182.9  m)  to  1.5  miles  (2.4  km)  nominal  width  (Morris  in  preparation).  The 
wider  reaches  actually  fluctuate  between  1.5  and  3.0  miles  (4.8  km)  in  width, 
depending  upon  local  Hood  conditions.  A  vertically  averaged  model  will  cap¬ 
ture  momentum  losses  due  to  expansions  and  contractions  of  the  river  channel, 
which  affect  the  progression  of  tides  and  Hood  waves.  Furlhcrmoro,  lateral 
momentum  transport  occurring  at  junctions  of  the  main  stem  with  tributaries 
will  be  well  represented.  These  characteristics  arc  important  for  the  ultimate 
computation  of  sediment  transport  and  water  quality,  especially  when  the  sedi¬ 
ment  load  or  water  chemistry  of  the  tributary  is  different  from  the  main  stem. 


Potential  Dimensionality 

The  LSJRE  has  been  categorized  as  a  highly  stratified  to  moderately 
stratified  estuary  by  Morris  (in  preparation),  who  reports  a  salt  wedge  of  vari¬ 
able  extent  in  the  LSJRE.  Vertical  stratification  due  to  salinity  or  temperature 
gradients  results  ir  vertical  velocity  gradients  and  secondary  circulation,  W'hich 
can  be  described  only  by  a  3-D  model  of  estuarine  hydrodynamics.  Secondary 
currents  can  also  result  from  How  around  channel  bends  and  from  wind  stress. 

A  3-D  mode!  represents  a  substantial  increase  in  sophistication  over  the 
2-D,  vertically  averaged  model.  Certain  modeling  objectives,  such  as  waste 
load  allocation  or  Hood  routing,  do  not  warrant  3-D  analysis.  However,  if  the 
model  results  arc  to  be  used  lor  modeling  salinity  stratification  and  the 
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response  of  phytoplankton  to  changing  salinity,  then  a  3-D  model  will  be 
required.  A  prudent  approach  would  be  to  select  a  2-D,  vertically  averaged 
hydrodynamic  model  that  has  a  3-D  counterpart  either  in  existence  or  actively 
under  development,  Since  most  3-D  models  are  built  upon  a  2-D  preliminary 
grid,  the  preparatory  work  of  grid  development,  boundary  condition  develop¬ 
ment,  and  data  gathering  will  carry  over  to  three  dimensions  when  the  need  for 
a  3-D  study  arises. 

Geometry 

Because  the  St.  Johns  River  is  bathymetrically  complex,  it  is  important  to 
define  the  complexities  accurately.  The  scales  of  flow  phenomena  vary 
because  of  the  dramatic  changes  in  width  and  depth  along  the  length  of  the 
river.  Either  a  stretched,  curvilinear,  boundary-lilted  finite  difference  grid  or 
an  unstructured  finite  element  grid  (or  finite  volume  grid)  would  be 
appropriate.  A  capability  for  rapid  local  refinement  of  the  grid  will  be 
advantageous  as  new  study  areas  are  defined. 

Wetting  and  Drying  Option 

Portions  of  the  St.  Johns  River  and  its  tributaries  are  adjacent  to  tidal  and 
freshwater  marshes.  Flooding  and  drying  of  marsh  areas  can  have  a  significant 
impact  on  the  volume  of  the  tidal  prism,  the  phasing  of  tidal  fluctuations,  and 
the  storage  volume  of  the  estuary.  As  part  of  the  LSJRE  ecosystem,  marshes 
may  play  an  important  role  in  filtering  sediments  and  water  quality  and  pro¬ 
vide  habitat  for  fish  and  wildlife.  A  model  of  LSJRE  hydrodynamics  should 
be  capable  of  handling  imcrmitlenl  flooding  of  marsh  environments  or 
representing  marsh  ‘'ora.  ;  schematically. 

Horizonte* '« M>*ViUient  Fluctuation  Terms 

Hydrodyn.  tic  circ  1  lion  in  the  LSJRE  is  controlled  largely  by  the 
complex  geomehy  ol  the  estuary.  Boundary  shear  is  important  around  expan¬ 
sions  and  contractions  of  the  river,  islands,  and  river  bends.  Horizontal 
turbulent  fluctuation  terms  are  thus  required  to  simulate  LSJRE  flows  ade¬ 
quately.  Wave  equation  models  that  ignore  horizontal  eddy  viscosity  are  not 
appropriate  lor  this  system. 

Vertical  Turbulent  Closure  Model 

An  advanced  turbulence  closure  model  such  as  the  Mcllor-  amadu  scheme 
is  recommended  for  the  3-D  hydrodynamic  model.  A  3-D  modeling  study  will 
be  usclii1  only  if  the  model  can  capture  the  stratification  behavior  of  the 
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estuary.  This  behavior  is  generally  controlled  by  the  vertical  turbulence  model. 
The  vertical  turbulence  closure  scheme  of  the  selected  model  should  be  well 
documented,  including  guidance  for  selecting  values  of  the  coefficients. 


Model  Maintenance 

The  models  described  in  this  report  are  actively  under  development.  The 
planning  period  for  development  of  the  SJKWMD  modeling  program  is  3- 
5  years,  ll  is  desirable  that  the  modeling  tools  developed  under  this  program 
be  useful  and  tip  to  dale  lor  several  years  beyond  the  development  program. 
Therefore,  it  is  prudent  to  choose  a  modeling  system  that  is  designed  to  take 
advantage  of  new  theoretical  developments  and  new  computer  technology. 
When  modeling  set  vices  arc  contracted,  provision  for  model  maintenance  and 
upgradability  should  be  negotiated. 

Given  the  above  ictpiircmcnL''.,  several  ol  the  models  reviewed  in  Chapter  6 
can  be  considered  most  appiopriate  for  the  LSJRL  system.  These  models  are 
(in  no  particular  order): 

a.  RMA-10  and  associated  models,  Ian  P.  King,  Wl-S, 

h.  Cl  I3D,  Y.  P.  Shong;  CH3D-WL.S,  Y,  1*.  Sheng  and  Lt.  II  Johnson 

e.  liCOM,  Alan  F.  Ulumbcrg  and  George  L  Mellor 

liach  of  these  models  lias  strengths  and  weaknesses  foi  the  LSJRLi  application. 
For  example,  model  a  is  the  only  one  of  the  three  that  lias  a  fully  implemented 
welling  and  drying  scheme.  Models  h  and  c  can  handle  marshes  schematically 
by  adding  "situ age  nodes."  Model  u  is  not  currently  linked  to  a  water  quality 
model  that  includes  phytoplankton.  Models  b  and  c  aie  each  linked  to  an 
advanced  walci  quality  model. 

Models  a  and  h  ate  public  domain  models  that  have  been  applied  to 
various  estuarine  flow  problems  by  the  authors  at  their  lespeclive  institutions, 
WFS  uses  both  of  these  models,  the  choice  of  which  is  made  according  to  the 
requirements  of  the  specific  application.  Model  c  is  proprietary  and  is 
available  by  contract  with  the  authors. 


Ctwpler  7  Recommended  Oiterla  foi  Selecting  u  HydiodynnitiicrSulinlty 


8  Hydrodynamic  Monitoring 
Requirements  to  Support 
Model  of  LSJRE 


Application  of  a  hydrodynamic/salinity  model  to  the  LSJRE  will  require 
field  data  to  set  boundary  conditions  and  to  perform  model  verification,  The 
minimum  data  requirements  arc  summarized  in  Table  1. 


Table  1 

Minimum  Data  Requirements 

Date  Category 

Date  Typo 

Sampling 

Locations 

Sampling 

Froquoncy 

... 

Period  ot 

Record 

Boundary 
conditions  and 
Initial  conditions 
(Continued) 

Tide 

elevation 

two  or  three 
locations  at  the 
ocean  boundary 

2.0-1 5.0  minutes 

Continuous  over 
spring-neap  tidal 
sequence,  plus 
long-term  records 
from  permanent 
tide  stations 
(existing  or  to  be 
established) 

"Velocity  or 
discharge 

Tributaries, 

point 

discharges,  and 

upstream 

boundary 

Daily 

Continuous 
record  ol  longth 
determined  by 
seasonality  of 
tributary  flows 

Salinity 

Minimally  0-  to 
12-point 
mousurements 
throughout  the 
estuary  (verti¬ 
cally  resolved) 

Duily 

Short-term 
records  (lor  initial 
conditions) 

T  ribu  lories, 
point 

discharges,  and 

upstream 

boundary 

Daily 

Continuous 
record  of  length 
determined  by 
sensonaity  of 
tributary  flow 

1  v  or  three 

locations  al 
ocean  boundary 
(vertically 
resolved) 

O.b  •  1  .0  hour 

Continuous  over 
one  spring-neap 
tidal  sequence 
for  each  season 

(Continued) 
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Table  1  (Concluded) 

Data  Category 

Data  Type 

Sampling 

Location* 

Sampling 

Frequency 

Period  of  Re¬ 
cord 

I  boundary 
conditions  and 
Initial  conditions 

Temperature 

Minimally  8-  to 
12-point 
measurements 
throughout  the 
estuary  (verti¬ 
cally  resolved) 

Doily 

Short-term 
records  (for  initial 
conditions) 

tributaries, 

point 

discharges,  and 

upstream 

boundary 

Daily 

Continuous 
record  of  length 
determined  by 
seasonality  of 
tributary  flows 

1  wo  or  three 
locations  at 
ocean  boundary 
(vertically 
resolved) 

0.5  -  1 .0  hour 

Continuous  over 
spring-neap  tidal 
sequence 

Wind  speed 
and  direction 

Point 

measurements 
throughout  the 
estuary 

Hourly 

Continuous 
record  of  length 
determined  by 
seasonality  ol 
local  meteorology 

Model 

verification  (two 
or  more 
independent 
data  sets 
required) 

Tidal 

elevation 

Minimally  6-  to 
12-point 
measurements 
throughout  the 
estuary  and  nt 
ocean  boundary 

“sinro 

minutes 

Continuous  over 
spring-neap  tidal 
sequence,  plus 
long-term  records 
from  existing 
permanent  tide 
stations 

Velocity 

Minimally  6-  to 
12-point 
measurements 
throughout  the 
estuary  and  at 
ocean  boundary 
(vertically 
resolved) 

~2.0~15.'0 

minutes 

Continuous  over 
spring -neap  tidal 
sequence 

Salinity 

Minimally  8-  to 

1 2-point 
measurements 
throughout  the 
estuary  and  at 
ocean  boundary 
(vertically 
resolved) 

0.5  •  1 .0  hour 

Continuous  over 
spring-neap  tidal 
sequence 

"Temperature 

Minimally  6-  to 

1 2-point 
measurements 
throughout  the 
estuary  and  at 
ocean  boundary 
(vertically 
resolved) 

0.5  -  1 .0  hour 

Continuous  over 
spring-neap  tidal 
sequence 
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9  Recommended  Approach 
for  LSJRE  Model 
Development  Program 


The  tasks  required  for  development  of  a  comprehensive  modeling  and 
monitoring  system  for  the  LSJRE  and  an  estimated  timeline  for  accomplishing 
these  tasks  are  given  in  Tables  2-11.  The  timeline  estimates  are  hypothetical, 
dependent  upon  the  satisfactory  completion  of  preliminary  tasks.  Furthermore, 
the  estimates  for  water  quality  and  sediment  modeling  system  development  are 
subject  to  modification  depending  upon  the  development  status  of  the  selected 
model. 
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Hydrodynamic  Model  Development  Program 


lescription 


Task 

H2 

Select  hydrodynamic  model  and  contract  tor  license,  model  maintenance,  and  consulting  services. 

H3 

Develop  limited  objectives  for  initial.  2-D  vertically  averaged  hydrodynamic  model  application  (time 
and  space  scales  to  be  simulated,  hydrologic  conditions  for  verification,  questions  to  be 
answered). 

H4 

Develop  2-D  vertically  averaged  hydrodynamic  grid  based  on  existing  maps  and  surveys. 

H5 

Obtain  preliminary  boundary  condition  data  from  historic  record. 

H6 

Perform  field  reconnaissance  to  verify  map  details. 

H7 

Perform  preliminary  simulations. 

H8 

Compare  model  to  historical  data. 

H9 

Analyze  and  document  preliminary  results. 

H10 

Determine  data  deficiencies  and  inadequacies.  Recommend  additions/deletions  In  hydrodynamic 
monitoring  network. 

H1 1 

Perform  further  model  simulations.  Continue  to  test  against  new  data. 

H12 

Analyze  and  document  secondary  results. 

Hi  3 

Develop  plans  for  future  2-D  vertically  averaged  model  applications.  Determine  if  existing  grid 
resolution,  validation,  and  analysis  system  is  adequate  tor  iuture  applications.  Develop  program 
of  model  maintenance  and  systematic  application  to  support  SJRWMD  management  objectives. 

H14 

Develop  limited  objectives  for  initial  3-D  hydrodynamic  model  application  (time  and  space  scales 
to  be  simulated,  hydrologic  conditions  for  verification,  questions  lo  bo  answered). 

HI  5 

Develop  3-D  hydrodynamic  grid. 

H16 

Develop  3-D  boundary  condition  data. 

H17 

Perform  preliminary  simulations 

HI  8 

Compare  model  to  historical  data. 

HI  9 

Analyze  and  document  preliminary  results. 

H20 

Determine  data  deficiencies  and  inadequacies  Recommend  additionsfdeletiuns  ill  hydrodynamic 
monitoring  network. 

H21 

Perform  further  model  simulations.  Continue  to  test  against  new  data. 

H22 

Analyze  and  document  secondary  results 

H23 

Develop  plans  for  future  3-D  model  applications.  Determine  if  existing  grid  resolution,  verification, 
and  analysis  system  is  adequate  for  future  applications  Develop  program  of  model  maintenance, 
enhancement,  and  systematic  application  to  support  SJRWMD  management  objectives.  j 

Develop  detailed  objectives  and  management  goals  lor  hydrodynamic  mooei  aovoiopmeru 
program. 

H2  Select  hydrodynamic  model  and  contract  for  license,  model  maintenance,  and  consulting  services. 

H3  Develop  limited  objectives  for  initial.  2-D  vertically  averaged  hydrodynamic  model  application  (time 

and  space  scales  to  be  simulated,  hydrologic  conditions  for  verification,  questions  to  be 
answered). 

H4  Develop  2-D  vertically  averaged  hydrodynamic  grid  based  on  existing  maps  and  surveys. 

H5  Obtain  preliminary  boundary  condition  data  from  historic  record. 

H6  Perform  field  reconnaissance  to  verify  map  details. 

H7  Perform  preliminary  simulations. 

H8  Compare  model  to  historical  data. 

H9  Analyze  and  document  preliminary  results. 

H10  Determine  data  deficiencies  and  inadequacies.  Recommend  additions/deletions  In  hydrodynamic 

monitoring  network. 

H11  Perform  further  model  simulations.  Continue  to  test  against  new  data. 

HI  2  Analyze  and  document  secondary  results. 

Hi  3  Develop  plans  for  future  2-D  vertically  averaged  model  applications.  Determine  if  existing  grid 

resolution,  validation,  and  analysis  system  is  adequate  tor  iuture  applications.  Develop  program 
of  model  maintenance  and  systematic  application  to  support  SJRWMD  management  objectives. 

HI 4  Develop  limited  objectives  for  initial  3-D  hydrodynamic  model  application  (time  and  space  scales 

to  be  simulated,  hydrologic  conditions  for  verification,  questions  lo  be  answered). 

HI  5  Develop  3-D  hydrodynamic  grid. 

HI 6  Develop  3-D  boundary  condition  data. 

H17  Perform  preliminary  simulations 

HI  8  Compare  model  to  historical  data. 

HI  9  Analyze  and  document  preliminary  results. 

H20  Determine  data  deficiencies  and  inadequacies  Recommend  additionsfdeletiuns  ill  hydrodynamic 

monitoring  network, 

H21  Perform  further  model  simulations.  Continue  to  test  against  new  data. 

H22  Analyze  and  document  secondary  results 

H23  Develop  plans  for  future  3-D  model  applications.  Determine  if  existing  grid  resolution,  verification, 

and  analysis  system  is  adequate  for  future  applications  Develop  program  of  model  maintenance, 
enhancement,  and  systematic  application  to  support  SJRWMD  management  objectives. 
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Tame  4 

Sediment  Transport  Model  Development  Program 


••cripUon 

veio^5«l?f#trcS]flc!?vnTrHrmSnaB«m8n^o3!^o^53!m6nnf«ns[poTmo3erSIvSIopifTM»ril^r3onwr,— 

52  Review  sediment  modeling  llte.ature. 

53  Select  sediment  modal  and  contract  for  license,  modal  maintenance,  modal  Interface  development  (if 
required),  and  consulting  tarvlcat. 

54  Uavelop  llmltad  ob|actlvas  for  Initial  sadlmant  modal  application  (time  and  space  scales  to  be  simulated, 
condlliono  for  verification,  questions  to  bo  answered). 

55  Initiate/expand  sediment  monitoring  program. 

SO  Obtain  sediment  chnierfe rlzatlnn  a>'d  boundary  condition  data. 

S7  Poilurm  preliminary  simulations. 

SO  Compare  model  results  to  prototype  to  the  extant  possible  witli  twisting  data. 

SO  Analyte  and  document  preliminary  result* 

510  Determine  dale  deficiencies  and  Inadsquacres.  Recommend  ur.ldlllt.<i>s/d«lutioris  In  sediment  monitoring 
network 

511  Perform  furlJroi  modol  simulations.  Contiriue  to  test  ayumM  now  data. 

51 2  Analyia  anti  document  secondary  results. 

513  Develop  pinna  foi  lulure  uidlmtmt  dynamic  model  applications.  Determine  II  oxlallng  gild  resolullon,  veilli- 
cation,  and  analyslA  system  Is  adequulo  lot  fuluri'  applications.  Develop  program  o!  model  mulnlenanoe 
arid  Byslanrallc  applicallon  lo  support  S.IRWMD  munageniant  objat.ilve*. 

. 1 - IT  I -TTTTnf  1IIIMBI  ,1  Ml  I  m 
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TablaS 

Tima  Una  for  Sadlmant  Transport  Modal  Dayajopmant  Program 


Y*ar  1 

Yaw  2 

Yaw  3 

Yaw* 

Years 

sr 

S2 

—  > 

S3 

--> 

S4 

- > 

GS 

— 

S6 

SY 

SO 

-> 

SO 

-> 

S10 

-> 

S11 

sin 

- > 

S13 
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Ctiapter  9  Racommencied  Approach  for  LSJRt:  Modal  Davalopmont  Program 


Water  Quality  Model  Development  Program 

hri*  “  '  "Description  " 

BT" - 


™ ^^ev^pTfe’tailedobjedivVs^nd'niiBnagement  goals  tor  water  quality  mooei  development  program, 

02  Review  water  quality  modeling  literature. 

03  Select  water  quality  model  and  contract  for  license,  model  maintenance,  model  Interface  development 

(If  required),  and  consulting  services. 

04  Develop  limited  objectives  tor  Initial  water  quality  model  application  (time  and  space  scales  to  be 

simulated,  conditions  for  verification,  questions  to  be  answered). 

Q5  Inltlate/expund  water  quality  monitoring  program. 

QG  Obtain  water  quality  boundary  condition  data. 

Q7  Perform  preliminary  simulations. 

QU  Compare  model  to  the  extant  possible  with  existing  data. 

09  Analyze  and  document  preliminary  results. 

Q10  Determine  data  deficiencies  and  Inadequacies.  Recommend  ndditions/deletlons  In  water  quality 

monitoring  network 

Oil  Perform  furthet  model  simulations.  Continue  to  test  against  new  data. 

Q12  Analyze  and  document  secondary  results. 

Q13  Develop  plans  lor  future  water  quality  mood  applications.  Determine  It  existing  grid  resolution, 

verification,  and  analysis  system  is  adequate  tor  future  applications.  Develop  program  ot  model 
maintenance  arid  systumatlc  application  to  support  bJRWMD  management  objectives. 
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Table  7 

Time  Line  for  Water  Quality  Model  Development  Program 

Teak 

Year  1 

Year  2 

Year  3 

Year  4 

Year  6 

Q1 

02 

03 

04 

_ > 

- > 

-> 

Qo 

06 

07 

06 

— 

— 

- > 

Q11 

U 
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Hydrodynamic  Monitoring  Program  Development 


escripuon 


ermine  gauge  locations,  sampling  intervals,  maintenance  schedules  tor  nydrodynamtc  data  stations 
in  consultation  with  modelers  and  end  users. 

Develop  database  for  hydrodynamic,  sediment,  and  water  quality  field  data  analysis,  storage,  and 
.etrieval. 

Install  gauges. 

Collect  and  malyze  data. 

Review  p.ocedures;  coordinate  with  modeling  team;  revise  program  as  needed. 


Sediment  Monitoring  Program  Development 


dft-cripuon 


termine  parameters,  procedures,  gauge  locations,  sampling  intervals,  maintenance  schedules  for 
sediment  rintn  stations. 

SM2  Expand  and  maintain  database  for  hydrodynamic  and  water  quality  field  data  analysis,  storage,  and 

reticivnl 

SM3  Obtain  samplers  and  equipment  if  required. 

SM4  c  ‘oiled  and  analyze  data. 

SM5  Review  procedures;  coordinate  with  modeling  team;  revise  program  as  needed. 


Water  Quality  Monitoring  Program  Development 


ascription 


etermine  parameters,  procedures,  gauge  locations,  sampling  intervals,  maintenance  scFiedulesTor 
water  quality  data  stations. 

I  xpand  and  maintain  database  for  hydrodynamic  and  water  quality  field  data  analysis,  storage,  and 
retrieval. 

Obtain  samplers  and  equipment  if  required. 

Collect  and  analyze  data. 

Review  procedures;  coordinate  with  modeling  team;  revise  program  os  needed. 
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Time  Line  for  Development  of  Monitoring  Programs 


10  Other  Considerations 


Computer  Hardware  Requirements 

Two-dimensional,  vertically  averaged  models  for  estuarine  hydrodynamics 
are  currently  being  run  on  minicomputers  and  midsize  workstations.  The  costs 
of  epu  and  memory  in  this  class  of  machines  have  steadily  decreased  over  the 
past  several  years,  and  this  trend  is  expected  to  continue.  Turn-around  times 
for  workstation  batch  jobs  rival  the  speeds  of  shared-resource  computers 
because  the  new  chips  arc  extremely  fast  and  the  user  has  local  control  over 
job  processing.  The  advantages  of  a  workstation  environment  include: 

a.  Graphical  user  interlace. 

b.  Moderate  initial  cost. 

c.  Personal  control  of  job  schedules. 

cl.  Multitasking. 

Three-dimensional  models  arc  computationally  very  expensive  and  time- 
consuming.  Current  run  times  for  fully  nonlinear  3-D  estuarine  models 
running  on  minicomputers  arc  typically  about  one-half  of  real  time,  depending 
on  the  size  of  the  grid  and  the  solution  technique.  Supercomputers  such  as  the 
CRAY  Y-MP  are  used  for  long-term  3-D  simulations.  However,  workstation 
epu  power  is  gaining  rapidly  on  the  supercomputer  market.  Low-end,  multi¬ 
processor  machines  arc  also  being  developed  that  can  take  advantage  of 
vcctorizablc  codes.  With  these  advancements,  a  locally  owned,  powerful  work 
station  with  periodic  upgrades  is  likely  to  be  the  most  cost-efficient  solution 
for  estuarine  modeling.  At  the  current  rale  of  workstation  development,  an 
evaluation  of  computer  resources  and  equipment  upgrade  is  recommended 
every  other  year. 

In  addition  to  the  epu,  the  following  equipment  and  software  arc  required: 

a.  Digitizer. 

b.  Graphics  monitor. 
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c.  FORTRAN  Compiler. 

d.  C  Compiler. 

e.  Plotting  software. 

f  High-quality  printer  lor  graphics. 

g.  Mass  data  storage  device. 

Graphical  User  Interface 

Estuarine  models  require  a  great  deal  of  spatial-  and  time-varying  data. 

The  models  generate  3-D  time-series  of  the  hydrodynamic,  sediment  transport, 
and  water  quality  parameters.  A  graphical  user  interface  is  necessary  to 
organize  information,  to  analyze  results,  and  to  gain  a  physical  understanding 
of  the  processes  described  by  these  data.  The  graphical  interface  should  be 
designed  to  manage  model  inpul/output,  visualize  3-D  data,  and  facilitate 
model  development  and  applications.  Severnl  graphical  user  interlaces  have 
been  developed  specifically  for  handling  data  sets  arising  from  geophysical 
flow  models.  A  graphical  user  interlace  may  be  purchased  from  a  third-party 
vendor  or  as  part  of  the  model  license.  As  with  other  components  of  the 
modeling  software,  a  contract  for  software  maintenance  and  updates  is 
recommended. 


Costs 

The  cost  of  developing  a  comprehensive  estuarine  flow  model  is  difficult  to 
estimate.  Costs  are  proportion.'))  to  the  degree  of  .sophistication  required  to 
answer  specific  management  questions.  Table  12  illustrates  "ballpark"  figures 
for  a  small-scale  model  development  plan  using  existing  codes.  The  costs 
cited  do  not  include  development  of  new  codes,  interfaces,  or  modeling  tech¬ 
nologies  that  may  be  required  to  address  certain  complex  problems.  The 
projected  costs  are  based  on  a  timetable  of  model  development  similar  to  that 
described  in  Chapter  9.  More  elaborate  plans  will  increase  program  costs. 

Costs  for  implementing  the  hydrodynamic,  water  quality,  and  sediment 
monitoring  programs  are  not  included  in  these  estimates.  It  is  assumed  that 
field  programs  will  be  conducted  by  local  agencies,  and  therefore  must  be 
budgeted  according  to  local  labor  and  equipment  resources.  A  comprehensive 
modeling  program  simply  will  not  succeed  without  coordination  between  the 
monitoring  program  and  the  modeling  program.  The  monitoring  costs  cited  in 
Table  12  reflect  the  efforts  of  (he  modeling  team  to  coordinate  the  field  effort, 
and  special  reconnaissance  related  specifically  to  modeling  needs. 
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Appendix  A  Notation 


B(x,z)  Width  of  the  estuary 

c  Celerity  of  a  gravity  wave 


D. 


One-dimensional  diffusion  coefficient  for  salinity 


D„Dy 


Turbulent  diffusion  coefficients  for  vertically  averaged 
salinity  concentration  in  the  x-  and  y-directions  {L2IT) 


D3?DyDz  Turbulent  diffusion  coefficients  for  salinity  in  the  x-,  y-,  and 
z-dircctions  (L2/T) 

F(s,T)  An  empirically  derived  relationship  describing  the  variation 
of  water  density  with  salinity  and  temperature 

g  Acceleration  due  to  gravity  acting  in  the  -z  direction  {LIT2) 

h  Bathymetric  depth  relative  to  datum  (L) 

it  Total  water  depth 

p  Local  pressuic  {MILT2) 

s  Salinity  concentration  {MIL3) 

S  Vertically  averaged  salinity  concentration  {MIL3) 

S  Cross-section  averaged  salinity  concentration  {MIL3) 

t  Time  (7) 

T  Temperature 

T  Cross-section  averaged  temperature 


Appendix  A  Notation 


u,v,w  Time-averaged  mean  velocities  in  the  x-,  y-,  and  z-directions 
{UT) 

u',v',w'  Turbulent  fluctuation  velocities  in  the  x-,  y-,  and  z-directicns 
(LIT) 

V  Depth-averaged  horizontal  velocity  in  the  x-direction  (LIT) 

XJ  Cross-section  averaged  velocity  (LIT) 

V  Depth-averaged  horizontal  velocity  in  the  y-direction  (LIT) 
x,y,z  Right-handed  Cartesian  coordinate  directions  ( L ) 

Longitudinal  tractive  forces  per  unit  volume  ( MIL 2T2) 

At  Time-step  for  simulation 

Axmin  Minimum  node  spacing 

Summation  of  tractive  and  body  forces  (excluding  gravity) 
per  unit  volume  in  the  x-,  y-,  and  z-directions  (MIL2!32) 

Turbulent  coefficient  of  viscosity  in  the  i'h  plane  in  the  fh 
direction,  also  called  eddy  viscosity  (MILT) 

£  Free  surface  elevation  relative  to  fixed  datum  (usually  mean 
sea  level)  (L) 

Source  (+)  or  sink  (-)  for  salinity  (MIL3!) 

0^  Cross-section  sources  and  sinks  for  salinity  (M/L3T) 


H  Viscosity  of  water  (MILT) 

P  Density  ol  water  (MIL3) 

p0  Reference  density  of  water  (M/L3) 

Tjj  Reynolds  stress  component  on  the  i'h  plane  in  the  fh 

direction,  for  i  =  \,  2,  3  and  j  -  1,  2  corresponding  to  the  x-, 
y-,  and  z-directions  (MILT2) 

Y  Cross-section  averaged  Reynolds  stress  in  the  longitudinal 
direction 


A2 
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